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Summary 


The  -transition  boundaries  for  a  triatomic  gas  between  various 
types  of  shock-wave  reflection  configurations — regular  reflection  (RR), 
single  Mach  reflection  (SMR) ,  complex  Mach  reflection  (CMR) ,  and  double 
Mach  reflection  (DMR) — exhibited  by  a  plane  shock  impinging  on  a  compressive 
corner  in  a  shock  tube  were  studied  experimentally  using  CO2  as  the  test  gas. 
The  shock  Mach  numbers  ranged  from  1.8  to  10.2,  and  the  initial  pressure  from 
3.5  to  80  torr  at  room  temperature.  The  regions  corresponding  to  the  four 
types  of  reflection  were  obtained  in  terms  of  the  shock  Mach  number^  Ma  jand 
the  wedge  angle,(6w.  The  resultant  transition  map  in  the  Ma  0w  plane  will 
be  of  value  to  future  researchers  in  the  field  as  an  aid  in  designing  their 
experiments.  -  5  m 

Excellent  agreement  was  obtained  between  the  transition  boundary 
from  CMR  to  DMR  according  to  the  criterion  proposed  by  Ben-Dor  §  Glass  (7a,  b) 
that  the  flow  downstream  of  the  reflected  shock  should  be  sonic  relative  to  the 
kink  of  the  reflected  shock  for  the  transition  to  occur  and  the  results  of  the 
present  experiments.  This  fact  will  further  confirm  the  validity  of  their  cri¬ 
terion. 


The  density  contours  of  the  flow  fields  were  obtained  by  means 
of  infinite-fringe  interferograms,  some  of  which,  particularly  those  of 
CMR  and  DMR,  were  hitherto  unavailable  in  the  literature.  They  may  prove 
valuable  for  comparison  with  the  results  of  numerical  simulation  of  the 
phenomena,  which  is  currently  the  subject  of  intense  research  at  various 
institutions.  *. 

The  analysis  of  shock-polar  diagrams  for  the  present  experiments 
appear  to  cast  doubts  on  the  validity  of  the  hypothesis  made  by  Henderson 
6  Lozzi  (29)  that  an  apparent  RR  beyond  the  von  Neumann  criterion  was 
actually  a  nascent  DMR  which  was  too  small  to  be  identified  as  such. 

It  seems  reasonable  to  conclude  from  the  results  of  the  present 
experiments  and  of  Ben-Dor  $  Glass  (7a,  b)  that  a  perfect-gas  model  with  an 
appropriate  constant  y  should  be  used  in  determining  the  transition  angles 
between  various  types  of  shock  reflection  for  shock  Mach  numbers  at  least 
up  to  10,  which  is  the  upper  limit  of  the  present  experiments. 
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complex  Mach  reflection 
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1.  INTRODUCTION 


When  a  plane  shock  wave  impinges  on  a  sharp  wedge,  it  exhibits 
different  patterns  of  reflection,  depending  on  the  strength  of  the  shock 
and  the  angle  of  incidence.  The  reflection  patterns  may  be  broadly  divided 
into  two  types:  regular  reflection  (RR)  and  Mach  reflection  (MR).  The 
latter  may  further  be  divided  into  three  subvarieties:  single  Mach  re¬ 
flection  (SMR) ,  complex  Mach  reflection  (CMR) ,  and  double  Mach  reflection 
(DMR)  (see  Fig.  1). 

The  shock  transition  between  various  types  of  reflection  patterns 
has  been  studied  extensively  in  the  past  not  only  as  an  important  gasdynamic 
phenomenon  per  se,  aiming  at  understanding  the  nature  of  reflection,  refraction, 
and  diffraction  of  shock  waves,  but  also  for  its  implications  and  utilities 
for  military  purposes,  such  as  the  prediction  of  shock  loading  on  structures 
and  optimal  heights  of  burst  for  weapons  [see  e.g.,  Glass  (22)].  Although 
the  studies  have  predominantly  been,  and  still  are,  concerned  with  the 
transition  between  RR  and  MR,  the  interest  in  the  transition  criteria  among 
other  types  of  reflection  has  been  increasing  in  recent  times;  especially 
with  the  advent  of  sensitive  pressure  transducers,  some  insight  into  the 
flow  properties  behind  the  reflected  shock,  which  are  perhaps  of  the  most 
practical  significiance,  has  been  obtained  (Merritt  (45),  Bazhenova  et  al  (4), 
Bertrand  (13)]. 

The  objectives  of  the  present  studies  are: 

1)  To  delineate  the  regions  corresponding  to  various  types  of  reflection 
in  carbon  dioxide,  in  terms  of  the  shock  Mach  number  Ms  and  the  actual 
wedge  angle  6W. 

2)  To  investigate  the  nature  of  the  secondary  Mach  shock*  which  characterizes 
DMR. 

3)  To  study  the  density  distribution  of  the  flow  fieJd  by  means  of 
infinite-fringe  interferograms . 

The  present  study  is  a  direct  continuation  of  the  work  by  Ben-Dor 
§  Glass  (7a, b)  who  studied  the  transition  criteria  of  shock  reflection  in 
nitrogen  and  argon,  extending  the  work  of  Law  f,  Glass  (43). 

The  transition  boundary  maps  depicting  the  regions  of  RR  and  MR 
have  appeared  in  previous  works:  e.g.  Courant  5  Friedrichs  (19),  Bleakev  5 
Taub  (14),  Polacheck  f,  Seeger  (50),  Syschikova  fi  Krassovskaya  (60). 

Bazenova  et  al.  (3)  were  the  first  to  attempt  to  delineate  the  regions 
corresponding  to  other  types  of  reflection  than  those  of  RR  and  MR.  Kith 
the  aid  of  pressure  transducers  embedded  in  the  reflecting  surface,  they 
obtained  over  a  limited  range  of  shock  Mach  number  the  transition  boundary 
between  CMR  and  DMR  experimentally.  Ben-Dor  fi  Glass  (7a, b)  proposed  a  new 
criterion  for  the  transition  from  CMR  to  DMR  that  the  flow  behind  the  re¬ 
flected  shock  be  sonic  with  respect  to  the  kink  of  the  reflected  shock,  and 
delineated  the  regions  for  four  known  types  of  reflection  —  RR,  SMR,  CMR, 
and  DMR  for  perfect  and  equilibrium  nitrogen  and  argon  over  the  shock 


*See  Sec.  2.2.3  for  the  description  of  the  secondary  Mach  shock. 
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Mach  number  range  from  1.0  to  8.0 


The  0  -  P/P q  shock-polar  diagram  is  provided  for  every  photograph 
of  the  present  experiments  included  in  this  report  so  as  to  afford  an 
overall  perspective  of  the  phenomena  under  investigation,  which  is  not 
visible  in  the  photograph.  For  example,  the  shock  polars  corresponding 
to  the  incident  and  reflected  shocks,  or  the  I-  and  R-polars,  together  with 
the  Mach  numbers  of  the  oncoming  flows  to  the  respective  shocks,  visualize 
all  the  states  achievable  downstream  of  the  shocks  in  terms  of  the  two 
parameters,  the  pressure  ratio  across  the  shock  P/Pq  and  the  flow  deflection 
angle  ■'< .  Moreover,  one  of  more  intersections  of  the  I  and  R-polars  are  the 
possible  solutions  of  the  two-shock  and  three-shock  relations,  so  that  the 
actual  states  of  the  angular  regions  around  the  confluence  point  of  the 
shocks  seen  in  the  photograph  are  indeed  represented  by  one  of  the  inter¬ 
section  points. 

It  should  be  borne  in  mind  that  these  0  -  P/Pq  shock -polar 
diagrams  are  the  pictorial  representations  of  the  two-shock  and  three- 
shock  interactions  in  the  frame  of  reference  with  its  origin  attached  to 
the  confluence  point  of  the  shock,  in  which  the  shocks  are  thus  made 
stationary. 

2.  CLASSIFICATION  OF  TYPES  OF  SHOCK  REFLECTION 

As  our  major  objective  is  to  delineate  the  regions  corresponding 
to  various  types  of  reflection  configurations  in  pseudo-stationary  flow, 
in  terms  of  the  shock  Mach  number  and  the  actual  wedge  angle,  it  is  desirable 
to  have  precise  ideas  of  what  is  meant  by  regular,  single  Mach,  complex  Mach, 
and  double  Mach  reflections  (Fig.  1)  and  how  they  are  determined. 

2. 1  Regular  Reflection 


Referring  to  Fig.  2a,  when  the  wave  angle  4>  is  sufficiently  small, 
the  intersection  point  P  of  the  incident  shock  I  and  the  reflected  shock  R 
occurs  on  the  wall,  and  moves  along  the  wall  surface.  In  a  coordinate 
system  attached  to  P,  the  initial  flow  uQ  along  the  wall  is  deflected  toward 
the  wall  by  the  incident  shock,  and  the  deflected  flow  u^  generates  the 
reflected  shock  R  which  turns  the  flow  uj,  so  that  U2  is  again  parallel 
to  the  wall.  Except  for  one  particular  value  $0  of  the  wave  angle,  however, 
the  reflection  is  not  specular,  i.e.  the  angle  of  reflection  w  is  not  equal 
to  the  wave  angle  j>0.  (The  value  of  <t>0  is  39.2°  for  y  =  1.4,  and  40.3°  for 
Y  =  1.29.)  For  j>o  greater  than  $0.  w  is  greater  than  4>0;  and  for  40  smaller 
than  $0,  w  is  less  than  4>0  [see  Polachek  §  Seeger  (50)]. 

In  the  coordinate  system  attached  to  the  reflection  point  P  (the 
intersection  point  of  the  two  shocks  and  the  wall),  the  flow  configuration 
is  made  stationary,  so  that  the  simple  two-shock  theory  [von  Neumann  (61), 
Bleakney  6  Taub  (14),  Polachek  5  Seeger  (50)]  may  be  used  to  find  the  flow 
properties  of  the  angular  regions  around  P.  Since  the  velocity  with  which 
P  moves  along  the  wall  is  very  high  for  small  values  of  wave  angle  <J>0>  t'1e 
flow  U2  behind  the  reflected  shock  R  is  in  general  supersonic  relative  to 
P,  except  for  the  narrow  range  of  <j>0  between  the  detachment  angle  (the 


limiting  wave  angle  beyond  which  RR  is  not  possible  according  to  the 
detachment  criterion,  see  Sec.  3  below,  and  the  sonic  angle  (the  wave 
angle  at  which  112  is  just  sonic  relative  to  P) ,  over  which  U2  is  sub¬ 
sonic.*  The  sonic  angle  is  also  referred  to  as  the  catch-up  angle,  since 
the  sonic  pulse  from  the  corner  can  catch  up  to  P  at  this  wave  angle. 

As  can  be  seen  from  the  wave  angle-flow  deflection  angle  diagrams  in 
Fig.  A2  in  Appendix  A,  the  sonic  angle  is  3°  to  4.5°  smaller  than  the 
detachment  angle  for  Mq  ^  2,  and  the  difference  quickly  decreases  with 
the  increases  in  Mq. 

If  the  flow  U2  behind  the  reflected  shock  is  supersonic  with 
respect  to  the  reflection  point  P,  the  sonic  signal  from  the  corner  cannot 
reach  P,  so  that  there  will  be  a  region  behind  R  in  which  the  thermodynamic 
state  of  the  gas  is  uniform;  as  a  result,  the  segment  of  R  delineating  this 
region  is  straight. 

Figure  2b  is  the  flow-deflection  angle-pressure  ratio  shock  polar 
diagram  of  regular  reflection  as  depicted  in  Fig.  2a.  Sonic  points,  of 
the  I -polar  and  S2  of  the  R-polar,  represent  the  condition  that  the  flows 
behind  the  incident  shock  I  and  the  reflected  R,  respectively,  are  just 
sonic.  Thus,  a  sonic  point  divides  the  polar  into  two  branches,  one  lying 
below  it  (usually  referred  to  as  the  weak  branch),  along  which  the  flow 
behind  the  shock  is  supersonic,  and  the  other  lying  above  it  (usually  re¬ 
ferred  to  as  the  strong  branch),  along  which  the  flow  behind  the  shock  is 
suhsonic.  Note  that  the  R-polar  meets  the  P-axis  below  the  intersection 
point  a  of  the  strong  branch  of  the  I -polar  and  the  P-axis,  and  that  there 
are  consequently  only  two  possible  solutions,  denoted  by  points  2  and  2'  in 
Fig.  2b.  Although  a  rigorous  theoretical  proof  is  yet  to  be  given,  it  has 
amply  been  demonstrated  experimentally  in  the  past  that,  in  case  two  points 
equally  satisfy  the  boundary  condition  that  0j  +62  =  0,  where  0^  and 
are  the  flow  deflection  angles  in  the  regions  behind  the  incident  and 
reflected  shocks,  respectively,  it  is  always  the  one  on  the  weak  branch 
that  is  actually  observed.  With  solution  point  2'  thus  eliminated  as 
physically  unlikely,  point  2  is  the  only  possible  solution.  Thus,  physical 
state  of  the  region  behind  the  reflected  shock  seen  in  the  shadowgraph  of  a 
typical  RR  in  Fig.  3a  is  represented  by  solution  point  2  in  Fig.  3b. 

However,  as  will  be  explained  in  detail  in  Sec.  3  below,  there 
are  two  conflicting  criteria--usually  referred  to  as  the  von  Neumann  and 
detachment  criteria--concerning  the  transition  from  RR  to  MR.  As  a  result, 
classification  of  the  reflection  type  becomes  subject  to  debate  in  the  case 
where  the  R-polar  meets  the  P-axis  above  the  intersection  a  of  the  I-polar 
and  the  P-axis,  see  Figs.  4b,  5b,  and  6b;  for  then,  excluding  solution 
point  2'  which  is  not  physically  meaningful,  there  now  exist  two  possible 
solutions  denoted  by  2N  (N  for  Neumann  criterion)  and  2p  (D  for  detachment 
criterion).  According  to  the  von  Neumann  criterion,  the  reflections  in 
Figs.  4a,  5a,  and  6a  are  Mach  reflection,  while  according  to  the  detachment 
criterion,  they  are  still  regular  reflection,  so  that  the  states  of  the 
regions  behind  the  reflected  shocks  (and  the  Mach  stems  in  the  case  of  the 
von  Neumann  criterion)  are  represented  by  2N  and  2n  respectively,  in 
Figs.  4b,  5b,  and  6b.  '  L 

To  the  accuracy  of  the  resolution  (0.2mm)  of  the  interferograms 
*See,  e.g..  Fig.  6  of  Bleakney  and  Taub  (14). 
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in  the  present  experiments,  reflections  in  Figs.  4a,  5a,  and  6a  are  regular. 
(For  further  discussion,  see  Sec.  3.) 

2.2  Mach  Reflection 


As  the  wave  angle  <j>Q  is  increased,  the  intersection  point  of  the 
incident  and  reflected  shock  occurs  off  the  wall,  and  is  connected  to  the 
wall  by  a  third  shock  M  (see  Fig.  7).  The  resultant  irregular  reflection 
configuration  is  called  Mach  reflection,  and  the  third  shock  is  referred  to 
as  the  Mach  stem. 

Since  the  confluence  consisting  entirely  of  three  shocks  is 
impossible  in  a  two-dimensional  flow  [Courant  5  Friedrichs  (19),  Breed  (17)], 
there  is  always  a  fourth  discontinuity  S  in  the  flow,  'across  which  entropy, 
density,  temperature,  and  flow  velocity  are  discontinuous,  bpt  pressure  and 
the  flow  direction  are  continuous,  and  is  referred  to  as  the  slipstream, 
contact  discontinuity,  or  vortex  sheet.  The  confluence  point  T  of  the  four 
discontinuities  is  called  the  triple  point. 

2.2.1  Single  Mach  Reflection 

This  configuration  occurs  when  the.  flow  U2  behind  the  reflected 
shock  R  is  subsonic  relative  to  the  triple  point  T,  and  is  identified  by 
the  curved  reflected  shock  over  its  entire  length  up  to  T  (see  Fig.  8a) . 

This  configuration  is  always  observed  when  the  incidence  is 
near  glancing  (<(>o  =  90°).  In  general,  in  the  case  of  SMR,  the  flow  angles 
to  the  incident  shock  and  the  Mach  stem  differ  only  slightly,  and  hence 
the  reflected  shock  is  weak.  As  a  result,  it  is  difficult  to  pinpoint  the 
triple  point  T  exactly  and  to  determine  the  shock  angles  at  T. 

In  most  cases,  the  shock  angles  at  T  are  measured  by  a  conventional 
protractor  or  some  modification  thereof;  however,  it  is  interesting  to  note 
that  Kawamura  §  Saito  (40)  determined  the  slopes  of  the  shocks  by  differenti¬ 
ating  the  5th-order  polynomials  approximating  the  traced  images  of  the  shocks. 

Figure  8b  shows  an  example  of  the  shock  polars  for  SMR,  which  is 
characterized  by  the  intersection  point  of  the  I-  and  R-polars  occuring  on 
the  subsonic  branch  of  the  R-polar.  Thus,  the  flow  u2  behind  the  reflected 
shock  R  is  subsonic  relative  to  the  triple  point  T.  Figure  9a  shows  the 
finite-fringe  interferogram  of  SMR.  In  general  the  flow  uj  behind  the 
incident  shock  I  is  only  slightly  above  sonic  in  SMR,  so  that  the  reflected 
shock  tends  to  be  nearly  normal  to  u^.  In  the  shock  polar  diagram,  this 
fact  is  manifested  by  the  smallness  of  the  R-polar  in  Fig.  9b,  for  which 
Mj  =  1.05.  Note  the  minute  R-polar  near  the  sonic  point  Sj  of  the  I-polar. 

At  the  limiting  flow  angle,  or  the  "extreme  sonic  angle" 

[Smith  (56)],  the  R-polar  degenerates  to  a  point  at  the  sonic  point  S,  of 
the  I-polar;  and  a  reflected  shock  is  no  longer  possible  beyond  this  limit, 
inasmuch  as  no  stationary  shock  can  exist  in  subsonic  flow.  This  means  that, 
in  the  polar  diagram,  no  R-polar  can  be  drawn  in  the  strong  branch  of  the 
I-polar  where  the  flow  is  subsonic. 
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Almost  laminar  isopycnics  along  the  reflected  shock  in  Fig.  10a 
show  in  the  most  striking  manner  that  the  expansion  wave  which  appears  in 
the  middle  of  the  reflected  shock  envelope  in  the  density  contours  obtained 
numerically  by  Schneyer  (51)  (see  Fig.  C-3c  in  Appendix  C)  is  indeed 
"spurious"  as  Schneyer  himself  remarked.  The  isopycnics  in  Fig.  10a  are 
evaluated  in  Fig.  10c,  using  the  frozen-gas  relations  with  y  =  1.290. 

Figure  11  illustrates  the  superposition  method  for  evaluating  band  dis¬ 
placement  for  the  case  comparative  to  Fig.  10a. 

2.2.2  Complex  Mach  Reflection 

When  the  flow  U2  behind  the  reflected  shock  R  is  supersonic 
relative  to  the  triple  point  T,  R  is  made  straight  at  T  [Kawamura  (39)]. 

The  resultant  Mach  reflection  configuration  (Fig.  12a)  with  a  straight 
reflected  shock  R  at  the  triple  point  T  is  referred  to  as  complex  Mach 
reflection  to  distinguish  it  from  single  Mach  reflection  which  has  a 
curved  reflected  shock  over  its  entire  length. 

It  should  be  borne  in  mind  that  there  is  nothing  intrinsically 
"complex"  about  this  type  of  reflection  in  particular,  as  the  term  may 
imply.  Because  of  the  straight  reflected  shock  in  CMR  and  DMR,  the  flow 
around  the  triple  point  is  uniform  and  indeed  simpler  in  CMR  and  DMR  than 
in  SMR  where  the  curvature  of  the  reflected  shock  R  inevitably  introduces 
non-uniformity  in  the  flow  behind  R. 


Figure  12b  is  a  typical  shock  polar  diagram  of  CMR,  which  is 
identified  by  the  intersection  of  the  I-  and  R-polars  occuring  in  the  super¬ 
sonic  branch  of  the  R-polar;  so  that  the  flow  behind  the  reflected  shock  is 
supersonic  relative  to  the  triple  point.  The  smaller  the  wedge  angle,  the 
closer  the  angle  of  incidence  to  the  extreme  sonic  angle,  and  hence  the 
closer  the  R-polar  to  the  sonic  point  S  of  the  I -polar. 

For  lower  range  of  Ms,  the  transition  from  the  straight  shock  R 
to  a  curved  bow  wave  may  take  place  smoothly  without  a  so-called  kink  K 
shown  in  Fig.  12a.  Figure  13a  illustrates  such  a  case.  As  Mg  is  increased, 
the  juncture  of  R  and  the  curved  bow  wave  becomes  more  and  more  pronounced 
as  can  be  seen  in  the  subsequent  photographs.  Note  the  Mach  waves  along 
the  wedge  surface  in  Fig.  13a,  which  are  generated  at  the  very  fine 
machining  striation  on  the  reflecting  surface.  The  density  of  each  isopycnic 
in  Fig.  14a  is  evaluated  in  Fig.  14c,  using  the  frozen-gas  relations  with 
y  =  1.290.  Figure  15a  is  a  shadowgraph  of  complex  Mach  reflection.  The  two 
plumb  lines  visible  in  this  and  other  photographs  are  1.0  cm  apart. 

The  dark  spots  visible  in  the  region  behind  the  curved  bow  wave,  one  large 
and  blurred,  and  the  other  small  and  distinct,  were  owing  to  superficial 
damages  on  the  lenses  of  the  interferometer  and  the  test-section  window 
glass.  Figure  15c  shows  the  Mach  numbers  of  the  flow  along  the  reflecting 
surface,  measured  by  means  of  the  Mach  waves  enamating  from  the  reflecting 
surface.  Considering  the  foot  of  the  Mach  stem  to  be  normal  to  the  wedge 
surface,  the  Mach  number  of  the  induced  flow  immediately  behind  it  was 
calculated  using  the  three-shock  theory,  and  the  predicted  value  1.88  is 
in  excellent  agreement  with  the  measured  value  1.9.  The  sudden  increase  in 
the  flow  Mach  number  downstream  of  the  Mach  stem  is  due  to  the  clockwise 


vorticity  of  the  flow  caused  by  the  interaction  of  the  flow  along  the 
wedge  surface  and  the  flow  along  the  slipstream  which  curls  up  near  the 
wedge  surface.  The  turbulent  nature  of  the  slipstream  is  clearly  visible 
in  Fig.  ISa.  As  shown  in  Figs.  16  and  17,  as  the  wedge  angle  decreases, 
the  angle  between  the  oncoming  flow  and  the  incident  shock  I  increases, 
causing  the  flow  velocity  behind  I  to  decrease.  Thus,  the  R-polars  in 
Figs.  16b  and  17b  are  located  close  to  the  sonic  points  S}  0f  the  I-polars. 
It  will  be  noticed  in  the  preceding  photographs  (Figs.  14a  through  17a)  that 
the  kink  K  becomes  more  and  more  pronounced  as  Mg  increases. 

2.2.3  Double  Mach  Reflection 


This  type  of  Mach  reflection  (see  Fig.  18a),  first  reported  by 
White  (64)  in  1952,  is  characterized  by  an  appearance  of  a  new  shock  M' , 
which  we  shall  call  the  secondary  Mach  shock*,  emanating  from  the  kink  of 
the  reflected  shock.  Consequently,  there  exist  two  systems  of  the  three- 
shock  interactions  in  the  flow,  hence  the  term  double  Mach  reflection. 

And  the  kink  now  becomes  the  secondary  triple  point  T'. 

For  lower  shock  Mach  numbers,  M'  extends  only  a  short  distance 
away  from  T',  as  can  be  seen  in  Fig.  19a  which  shows  a  nascent  secondary 
Mach  shock  M' .  As  well,  in  Fig.  19a,  the  derivative  of  the  slope  of  the 
shocks  at  T'  appears  to  be  still  continuous,  and  the  entire  straight 
segment  of  the  primary  slipstream  S  is  distinctly  visible  as  the  centre  of 
the  vorticity  of  the  flow  at  the  juncture  of  S  and  the  reflecting  surface 
is  some  distance  away  from  S,  but  its  influence  is  strong  enough  to  cause 
the  Mach  stem  to  slightly  bulge  upstream.  The  reflected  shock  R  from  the 
primary  triple  point  T  in  Fig.  19a  is  almost  parallel  to  the  trajectory 
path  of  T,  inclined  at  the  angle  x  to  the  reflecting  surface  (Fig.  18a). 

In  Fig.  20a,  weak  though  it  may  be,  a  spoon-shaped  M'  spans  between  T' 
and  the  primary  slipstream  S,  intersecting  the  latter  almost  perpendicularly. 
The  image  of  M'  at  T1  appears  as  though  M'  were  the  appendage  of  the  curved 
bow  wave,  no  change  in  the  shapes  of  the  two  shocks  being  discernible. 

A  sudden  change  of  the  shape  of  M'  near  S,  as  a  result  of  which  M'  has  a 
trowel-like  appearance,  may  be  construed  as  an  indication  of  complexity  of 
the  flow  downstream  of  the  reflected  shock  R.  The  image  of  the  primary 
slipstream  S  becomes  progressively  ill-defined  beyond  the  half-length, 
illustrating  the  fact  that  the  shear  flow  along  S  near  the  primary  triple 
point  T  turned  into  turbulence. 

Photographs  in  Figs.  21a  through  24a  were  taken  using  the  same 
wedge  of  27°  semi-vertex  angle,  and  are  in  the  increasing  order  of  Ms  so 
as  to  visualize  the  corresponding  gradual  growth  of  M' .  Figures  21a  shows 
an  incipient  M'.  (Dark  images  near  the  reflecting  surface  are  due  to  the 
smudges  on  the  test-section  window.)  Figure  22a  distinctly  shows  the  nature 


*This  terminology  was  chosen  so  that  each  shock  in  the  flow  field  may  be 
uniquely  referred  to,  even  though  it  differs  from  the  system  of  reference  adopted 
by  Ben-Dor  8  Glass  (7a, b)  who  considered  the  new  shock  under  discussion  as  the 
secondary  reflected  shock  Rj  by  drawing  analogy  between  the  gasdynamic 
discontinuities  I,  R,  M,  S  and  R,  R^,  Sj,  respectively,  see  Fig.  18c. 
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of  turbulence  along  S.  Because  M'  i«  very  weak,  the  secondary  slipstream 
is  entirely  invisible,  even  though  its  existence  is  guaranteed  theoretically 
as  previously  mentioned.  In  Fig.  23a,  M' ,  smooth  and  distinct,  connects  T' 
and  S,  meeting  S  somewhat  obliquely.  (Again,  the  smudges  on  the  window  are 
visible.)  To  the  accuracy  of  the  resolution  of  the  photograph,  the 
derivative  of  the  slope  of  the  shock  envelope  at  T'  is  continuous. 

With  further  increases  in  Ms,  M'  becomes  more  and  more  straight, 
and  its  intersection  point  with  S  shifts  downstream  towards  the  reflecting 
surface  as  evident  in  F.s,  24a.  Figures  25a  through  32a  show  cases  of 
strong  to  very  strong  incident  shocks.  In  Fig.  25a,  the  curvature  of  M' 
appears  to  be  concave  downstream,  while,  ir  the  previous  photographs,  it 
is  concave  upstream.  The  intersection  of  M*  with  S  occurs  nearly  at  the 
extreme  position  relative  to  S,  i.e.  near  the  reflecting  surface.  The 
angle  w  between  the  reflected  shock  R  and  the  trajectory  path  of  the 
triple  point  T,  or  the  angle  of  reflection  (see  Fig.  18a),  is  now  "negative" 
that  is  to  say,  R  is  located  below  the  trajectory  path.  Visible  in  the 
shadowgraph  of  Fig.  25a  is  the  secondary  slipstream  S'  enamating  from  T' . 
Note  the  crescent  image  in  the  lower  part  of  the  Mach  stem  M.  This  is  a 
manifestation  of  the  three-dimensional  character  of  the  bulge  of  M  as 
explained  below,  since  its  effect  is  more  apparent  in  the  later  photographs. 
As  in  Fig.  20a,  S  becomes  less  and  less  visible  along  its  length,  indicating 
the  gradual  transition  from  shear  flow  to  turbulence  along  S.  A  large  curl¬ 
ing  of  S  near  the  reflecting  surface  observable  in  the  previous  photographs 
has  drastically  diminished  its  size.  The  gray  furrow-like  image  along  the 
reflecting  surface  behind  the  curved  bow  wave  is  the  boundary  layer, 
enhanced  by  the  minute  striations--from  which  Mach  waves  are  generated-- 
on  the  surface.  Figure  26a  confirms  the  observation  made  earlier,  i.e. 

M*  is  strongest  at  T'.  This  is  borne  out  by  the  fact  that  the  isopycnics 
of  M'  become  denser  towards  I'.  A  high  degree  of  compression  of  the  gas 
near  the  foot  of  the  slipstream  Sis  made  manifest  by  the  cluster  of  dense 
isopycnics.  No  signs  of  the  secondary  slipstream  are  evident.  Note  the 
disturbance  immediately  ahead  of  the  curved  bow  wave.  This  disturbance, 
the  cause  of  which  is  yet  to  be  investigated,  appears  quite  distinctly  in 
some  interferograms ,  not  only  of  the  present  experiments  (Figs.  27a,  28a, 
and  29a  below,  for  example),  but  also  in  an  interferogram  in  CO2  included 
in  Law  6  Glass  (43)  who  used  the  same  apparatus  as  in  the  present  case. 

A  remarkable  fact  is  that  it  is  only  in  (Jib  that  the  disturbance  is  evident, 
and  not  in  any  other  types  of  gases.  The  density  contours  such  as  shown  in 
Fig.  27a  makes  one  wonder  whether  some  extraneous  disturbances  are 
influencing  the  flow  behind  the  curved  bow  wave.  Figure  27a  shows  the 
flow  downstream  of  the  Mach  stem  M  surging  up  from  the  reflecting  surface 
along  the  vortex  of  the  slipstream  toward  M.  This  is  the  flow  that  causes 
M  to  bulge  near  its  foot.  The  density  of  each  isopycnic  in  Fig.  27a  is 
evaluated  in  Fig.  27c,  using  the  frozen-gas  relations  with  y  =  1.290,  and 
the  results  obtained  by  the  finite-fringe  method  by  Ben-Dor  &  Glass  (7a, b) 
for  comparable  cases  in  N2  as  the  guide.  The  values  for  p/p0  assigned  to 
the  isopycnics  in  the  region  inside  the  curved  bow  wave  appear  to  be  rather 
large.  This  resulted  from  the  fact  that  the  variation  in  density  between 
two  consecutive  isopynics  in  an  infinite-fringe  interferogram  is  known 
(see  Appendix  C) ,  and  the  fact  that  the  values  for  p/p0  increase  steadily 
toward  downstream  according  to  the  isopynics  in  Ref.  7a, b.  Figures  28a 
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through  31a  are  of  the  experiments  using  the  identical  wedge  of  30°  semi- 
vertex  angle.  (The  Mach  waves  emanating  from  the  reflecting  surface  are 
particular  to  this  wedge.)  Due  to  the  constraints  on  the  experimental 
set-up,  the  initial  pressures  are  lower  for  large  Ms.  As  noted  above, 

Figs.  28a  and  29a  show  the  disturbance  immediately  ahead  of  the  curved 
bow  wave.  No  trace  of  the  secondary  slipstream  is  discernible  in  either 
figure.  In  Fig.  29a,  the  bulge  of  M  can  be  examined  in  detail.  The 
clockwise  flow  along  the  vortex  associated  with  the  slipstream  S  interacts 
with  the  induced  flow  behind  M;  as  a  result,  a  part  of  M  is  "pushed"  upstream, 
causing  a  bulge.  In  Fig.  29a,  the  bulge  is  seen  to  be  three-dimensional  in 
the  sense  that  the  extent  of  the  departure  from  the  planarity  of  the  front 
is  more  severe  in  the  region  near  the  centre-plane  of  the  shock  tube  than 
near  the  window,  as  a  close  examination  of  the  profile  of  M  will  show. 

(In  Fig.  25a  above,  one  sees  the  early  stage  of  this  effect.)  In  Figs. 

30a  and  31a,  the  concentric  rings  of  isopycnics  indicative  of  the  vorticity 
associated  with  the  slipstream  S  are  no  longer  visible;  the  slipstream 
appears  L-shaped,  instead.  The  number  of  isopycnics  visible  in  either 
interferogram  is  less  than  in  the  previous  ones,  due  to  the  decrease  in 
the  initial  density.  Figure  32a  shows  the  highest  Ms  achieved  in  the 
present  experiments  (Ms  =  io. 17) >  which  also  results  in  DMR.  The  bulge 
of  M  appears  two-dimensional  in  this  case.  Note  the  somewhat  peculiar 
shapes  of  isopycnics  downstream  of  M,  indicating  the  complexity  of  the  flow 
field.  The  density  of  each  isopycnic  in  Fig.  32a  is  evaluated  in  Fig.  32c, 
using  the  frozen-gas  relations  with  y  =  1.290. 

It  is  not  possible  to  distinguish  between  CMR  and  DMR  on  the 
basis  of  the  shock  polar  diagrams,  since  the  flow  behind  the  reflected  shock 
is  supersonic  in  either  case. 

3.  TRANSITION  FROM  RR  TO  MR 


As  mentioned  in  Sec.  2.1  above,  there  are  two  main  criteria 
for  the  transition  from  RR  to  MR.  One,  attributed  to  von  Neumann  (61), 
requires  that,  when  the  flow  angle  to  the  incident  shock  is  increased 
while  maintaining  the  shock  strength  constant,  the  changes  in  the  pressures 
in  the  regions  around  the  confluence  point  of  the  shocks  should  also  be 
gradual  at  transition  from  RR  to  MR.  In  other  words,  the  regions  will 
remain  in  "mechanical  equilibrium"  during  the  transition.  According  to 
this  criterion,  therefore,  the  pressure  in  the  region  immediately  behind 
the  reflected  shock  at  the  confluence  point  will  not  suffer  a  sudden  drop 
when  the  transition  occurs,  the  phenomenon  predicted  by  the  detachment 
criterion  discussed  below.  [This  is  most  easily  understood  by  examination 
of  the  flow  deflection  angle  -  pressure  ratio  (6  -  P/P^)  shock  polar 
diagrams,  see  Figs.  4b,  5b,  and  6b,  and  also  Fig.  2  of  Homung  et  al  (32).] 
The  required  condition  is  realizable  only  in  the  so-called  stationary  Mach 
reflection  configuration  (Fig.  33)  in  which  the  slipstream  is  parallel  to 
the  wall  and  the  Mach  stem  is  perpendicular  to  the  oncoming  flow  [von  Neumann 
(61),  Courant  §  Friedrichs  (19)].  Thus,  one  of  the  boundary  conditions  in 
this  configuration  is  that  the  flow  deflection  suffered  through  the  incident 
shock  be  exactly  compensated  by  the  re-deflection  through  the  reflected 
shock. (See  Refs.  19,  50,  61  for  further  discussion.) 
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The  other,  usually  referred  to  as  the  detachment  criterion, 
which  is  put  forward  by  standard  textbooks  on  gasdynamics  such  as 
Liepmann  5  Roshko  (44) ,  specifies  that  the  transition  from  RR  to  MR 
occurs  when  the  flow  deflection  angle  through  the  incident  shock 
exceeds  the  maximum  flow  deflection  angle  C>2  through  the  reflected  shock. 

The  flow  angle  <)>0  to  the  incident  shock  corresponding  to  the  von  Neumann 
criterion  is  less  than  that  corresponding  to  the  detachment  criterion  over 
the  entire  range  of  shock  strength,  except  at  one  critical  point  at  which 
the  two  values  coincide,  see  Fig.  34a.  (Note  that  the  effective  wedge 
angle  0^  in  Figs.  34a-e,  which  is  the  sum  of  the  actual  wedge  angle  6W 
and  the  trajectory  angle  x  of  the  triple  point,  is  the  complement  of  the 
flow  angle  <j>0  to  the  incident  shock,  i.e.  0^  =  tt/2 - <t>0) . 

In  addition  to  the  above  two  main  criteria,  a  third  one,  which 
may  be  termed  the  sonic  criterion,  has  been  suggested  [see,  e.g.  Kawamura  5 
Saito  (40) ] according  to  which  the  transition  from  RR  to  MR  is  expected  to 
occur  at  the  sonic  angle,  the  angle  of  incidence  at  which  the  flow  behind 
the  reflected  shock  is  just  sonic  relative  to  the  reflection  point.  The 
underlying  argument  is  that,  when  the  flow  behind  the  reflected  shock  is 
subsonic,  the  cor  ier  signal  ca/i  reach  the  reflection  point  P  causing 
disturbances  in  the  entire  flow  downstream  of  P,  and  the  transition  to 
MR  is  set  off  as  a  result,  since  by  assuming  the  MR  configuration,  the 
flow  is  relieved  of  the  high  pressure  behind  the  reflected  shock  and  hence 
is  more  "stable"  than  just  prior  to  the  transition  [see  e.g.  Jahn  (35)]. 

However,  owing  to  the  proximity  of  the  sonic  and  the  detachment 
angles,  it  is  not  possible  to  resolve  these  two  points  experimentally. 

In  the  region  in  the  angle  of  incidence-shock  strength  plane, 
or  the  4>0  _  £  plane  (or  equivalently  in  the  Ms  -  0^  plane,  as  done  in  this 
study) ,  overlapped  by  two  boundaries  corresponding  to  the  von  Neumann  and 
detachment  criteria,  either  type  of  reflection,  RR  or  MR,  is  theoretically 
possible.  So  far  as  strong  shocks  are  concerned,  however,  the  results  of 
shock-tube  experiments  by  Smith  (56),  Bleakney  §  Taub  (14),  Kawamura  6 
Saito  (40),  Henderson  Lozzi  (29),  Hornung  &  Kychakoff  (31),  Ben-Dor  §  Glass 
(7a, b),  and  others  have  given  evidence  in  support  of  the  detachment  criterion, 
while  the  experiments  by  Hendorson  §  Lozzi  (29),  Hornung  §  Kychakoff  (31), 
and  Hornung  et  al  (32)  appear  to  confirm  the  von  Neumann  criterion  in  the 
case  of  stationary  flow. 

In  light  of  this  apparent  inconsistency  of  the  transition  criteria 
between  the  two  types  of  flow,  Henderson  6  Lozzi  (29),  and  later, 

Syshchikova  5  Krassovskaya  (60) ,  postulated  that  the  transition  from  RR  to 
MR  occurred  at  the  boundary  corresponding  to  the  von  Neumann  criterion,  and 
that  the  apparent  RR  beyond  that  limit  was  actually  a  nascent  double  Mach 
reflection  that  was  too  small  to  be  identified. 

On  the  other  hand,  Hornung  et  al  (32)  suggested  an  explanation 
for  the  well-known  persistence  of  RR  beyond  the  limit  predicted  by  the 
detachment  criterion  in  terms  of  the  displacement  thickness  of  the  boundary 
layer  at  the  reflection  point  in  pseudo-stationary  flow.  (We  note  that 
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this  suggestion  is,  however,  not  without  difficulty;  for  as  mentioned  in 
their  paper,  previous  experiments  by  one  of  the  authors  had  shown  an 
increase  in  the  angle  of  reflection,  while  the  proposed  displacement 
effect  of  the  wall  should  have  resulted  in  a  decrease.) 

Thus,  there  exist  two  opposing  views  regarding  the  persisted 
regular  reflection:  according  to  Henderson  5  Lozzi  (29),  it  is  a 
nascent  double  Mach  reflection,  while  Hornung  et  al  (32)  claim  that  it 
is  still  a  regular  reflection. 

On  the  basis  of  the  appearance  in  the  photographs,  a  series  of 
present  experiments  done  by  using  a  wedge  of  50.4°  semi -vertex  angle 
(three  of  which  are  shown  in  Figs.  4,  5,  and  6)  have  been  classified  as 
regular  reflection.  According  to  the  hypothesis  of  Henderson  §  Lozzi  (29) , 
however,  they  are  all  double  Mach  reflection;  and  the  states  behind  the 
reflected  shock  and  Mach  stem  map  into  the  point  2^  in  the  accompanying 
shock  polar  diagrams.  (The  poii t  2q  corresponds  to  regular  reflection.) 

According  to  our  three  •:.!  ,  .;k  solution  for  y  =  1.290,  the  flow 
angles  4>3  to  the  Mach  stems  corresponding  to  the  point  2N  in  Figs.  5  and  6 
are  88.5°  and  87.5°,  respectively,  lienee,  if  one  assumes  that 

X  =  9J°-*S 

then  the  predicted  angles  x  of  the  triple  point  trajectory  will  be  1.5°  and 
2.5°,  respectively,  and  therefore,  the  Mach  stem  should  definitely  be 
observable. 


The  two-shock  theory  gives  two  values  for  the  angle  of  reflection 
w  for  given  shock  strength  £  and  angle  of  incidence  <j>0;  the  higher  value  is 
termed  the  strong  solution  and  the  lower  one  the  weak  solution,  and  it  is 
the  latter  that  is  usually  observed  in  practice.  The  two  solutions  coincide 
at  a  critical  angle  of  incidence,  or  the  extreme  angle,  beyond  which  no 
solutions  exist.  Polacheck  5  Seeger  (50)  for  example,  show  that  the  solution 
curves  of  the  two-shock  and  three-shock  theories  for  a  given  shock  strength 
£  in  the  <J>  -  w  plane  intersect  at  a  single  point.  The  particular  value  of 

£  for  which  this  intersection  point  coincides  with  the  extreme  angle  of  the 
two-shock  curve  is  the  critical  shock  strength  that  separates  the  two  classes 
of  shocks,  weak  and  strong.  [As  shown  by  Polacheck  £  Seeger  (50),  this 
occurs  at  £  =  0.433  and  q0  =  41.4°  for  y  =  1.4.J  Kawamura  §  Saito  (40)* 
found  that  in  the  case  of  strong  shocks,  the  transition  from  RR  to  MR  was 
accompanied  by  a  sudden  decrease  in  the  angle  of  reflection  w  while  the 
transition  was  smooth  and  the  angle  of  reflection  increased  continuously  in 
the  case  of  weak  shocks. 


*There  has  been  a  certain  amount  of  uncertainty  concerning  the  numerical  values 
of  the  parameters  cited  in  Kawamura  §  Saito' s  paper,  steming  mainly  from  an 
apparent  misprint  of  the  value  of  the  Mach  number  of  the  incident  flow.  Molder 
(48)  recently  recalculated  these  quantities  for  y  =  7/5  and  5/3,  maintaining 
seven  significant  digits. 
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Here  we  disagree  with  Henderson  §  Lozzi  (30)  who  interpreted 
Kawamura  5  Saito's  findings  to  be  the  "tojj  [the  von  Neumann  criterion]  was 
the  criterion  for  strong  shocks  and  the  tog  [the  detachment  criterion]  was 
the  one  for  weak  shocks".  As  evident  from  the  inspection  of  the  0o-u  planes 
or  the  q-P  polar  diagrams,  the  transition  from  RR  to  MR  according  to  the  von 
Neumann  criterion  must  be  smooth  without  a  discontinuous  change  in  the 
angle  of  reflection,  while  the  detachment  criterion  predicts  a  sudden  drop 
in  the  angle  of  the  reflection  at  the  transition.  Thus,  the  results  of 
strong-shock  experiments  by  Kawamura  5  Saito  are  reconcilable  only  with 
the  detachment  criterion. 


4.  TRANSITION  FROM  CMR  TO  DMR 

Law  §  Glass  (43)*  suggested  that  DMR  occured  when  the  flow  U2 
(see  Fig.  8)  behind  the  reflected  shock  R  was  supersonic  relative  to  the 
wedge  surface.  Gvozdeva  et  al  (26)  suggested  that  the  DMR  configuration 
was  brought  about  as  a  result  of  the  excitation  of  the  internal  degrees 
of  freedom  of  the  gas  molecules  behind  the  shocks.  Semenov  et  al  (53) 
proposed  that  the  inception  of  the  secondary  Mach  shock  M'  (Fig.  21a)  was 
due  to  the  rotational  flow  which  developed  in  the  region  between  the 
primary  Mach  stem  M  and  the  reflected  shock  R  and  caused  the  primary  slip¬ 
stream  S  to  curl  up  into  a  vortex.  Gvozdeva  et  al  (26)  and  Bazhenova  et  al 
(4)  measured  the  pressure  on  the  reflecting  surface  of  the  wedge  by  using 
piezo-gauges  and  found  that  their  results  concerning  the  inception  of  DMR 
did  not  agree  with  the  criterion  suggested  by  Law  §  Glass  (43).  They 
considered  that  the  flow  U2  behind  the  reflected  shock  R  to  be  supersonic 
with  respect  to  the  primary  triple  point  T  as  the  necessary  (but  not 
sufficient)  condition  for  DMR,  and  suggested  that  the  appearance  of  the 
second  increase  in  pressure  on  the  reflecting  surface  from  the  oscillo¬ 
graphic  trace  as  a  way  to  distinguish  DMR  from  other  types  of  Mach  reflec¬ 
tion.  We  remark,  however,  the  difference  in  sensitivity  between  the  optical 
method  and  the  pressure-transducer  method  of  data  acquisition.  Figure  19  is 
a  case  in  point;  while  a  nascent  secondary  Mach  shock  is  clearly  visible  in 
the  shadowgraph,  it  is  doubtful  if  it  will  make  its  presence  known  to  the 
pressure  transducer  on  the  reflecting  surface,  which  after  all  is  capable 
of  registering  only  the  integrated  pressure  over  the  area  of  its  pressure- 
sensitive  surface. 

Thus,  it  seems  natural  that  the  two  boundaries  determined 
experimentally  by  an  optical  method  and  a  pressure- transducer  method  should 
disagree.  Ben-Dor  Glass  (7a, b)  suggested  that  the  transition  from  CMR  to 
DMR  occurred  when  the  flow  U2  behind  the  reflected  shock  R  (see  Fig.  12a) 


*In  the  paper  cited,  Law  $  Glass  acknowledge  the  suggestion  made  to  them  by 
Dr.  L.G.  Gvozdeva  of  the  Institute  of  High  Temperature,  Moscow,  U.S.S.R.  as 
the  basis  of  their  criterion. 
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was  supersonic  with  respect  to  the  kink  K  of  R.  They  compared  the  predicted 
regions  of  various  types  of  reflection  in  the  Ms  -  0^  plane,  where  Mg  is 
shock  Mach  number  and  9^  the  "effective  wedge  angle"  (the  sum  of  the  wedge 
angle  6W  and  the  angle  x  of  the  triple-point  trajectory),  with  experiments 
using  argon  and  nitrogen  over  the  range  of  M  from  1.9  to  8.  The  results 
for  N2  are  shown  in  Fig.  53  of  Ben-Dor  (6)  and  Fig.  18  of  Ben-Dor  5  Glass 
(7a).  In  spite  of  an  apparent  agreement  of  data  points  with  the  predicted 
regions  corresponding  to  equilibrium  nitrogen  as  shown  in  the  two  figures, 
if  one  plots  the  data  for  the  experiments  in  N2  listed  in  Table  5  of 
Ben-Dor  (6)  onto  the  Mg  -  9^  plane  in  Fig.  23  of  the  same  report,  one  will 
notice  that  there  are  some  points  that  fall  into  the  regions  overlapped  by 
the  two  boundaries  corresponding  to  a  perfect  diatomic  gas  y  =  1.4  and  to 
equilibrium  nitrogen  (e.g.  Exp.  #33,  42,  43,  44,  45,  46  of  Ref.  6).  The 
reason  for  the  fact  that  the  experimental  points  now  agree  with  the  regions 
corresponding  to  a  perfect  gas,  instead  of  those  corresponding  to  equili¬ 
brium  gas  may  be  that  the  theory  is  adequate  but  the  equilibrium-gas  analysis 
is  irrelevant  for  the  purpose  of  delineating  the  regions  since  the  shock 
waves  remain  frozen  (Appendix  E) . 

Since  the  criterion  for  the  ae...  iti 01.  from  CMR  to  DMR  put 
forward  by  Ben-Dor  5  Glass  (7a, b)  is  in  terms  of  the  Mach  number  of  the 
flow  behind  the  reflected  shock  with  reference  to  the  kink  K  which  itself 
is  moving  in  the  laboratory  frame  of  reference,  it  is  necessary  to  find  the 
relative  velocity  u  of  the  kink  with  respect  to  some  reference  point,  say, 
the  triple  point  T.  Thus,  Ben-Dor  §  Glass  (7a, b)  assumed  as  Law  6  Glass  (43) 
did  that  the  kink  K  moved  with  the  same  horizontal  velocity  as  the  induced 
flow  behind  the  incident  shock  I,  so  that  the  ratio  l/l,  where  l  and  L  are 
horizontal  distances  between  T  and  K,  and  T  and  the  vertex  of  the  wedge, 
respectively,  is  equal  to  the  density  ratio  across  the  incident  shock 
irrespective  of  the  wedge  angle  9W.  This  assumption,  however,  was  proved 
to  be  invalid  by  Bazhenova  et  al  (4)  who  measured  the  ratio  l/l  in  their 
experiments  and  found  that  for  a  given  incident  shock  Mach  number  the  ratio 
decreased  with  the  increases  in  the  wedge  angle. 

It  is  relevant  here  to  remember  that  along  the  transition  boundary 
between  CMR  and  DMR  (see  Figs.  34  and  35),  the  ratio  l/l  obviously  vanishes 
because  the  triple  point  T  and  the  kink  K  coincide,  and  hence  l  =  0,  on  that 
boundary. 


In  the  present  experiments,  the  measurements  were  made  of  the 
ratio  l/l  for  three  incident  shock  Mach  numbers  5,  7,  and  8;  the  results 
are  shown  in  Fig.  36. 

In  general,  the  pressure  in  regions  4  and  5  behind  the  bow  wave 
(see  Fig,  18a)  is  greater  than  that  in  regions  2  and  3  behind  the  reflected 
shock  R,  as  evident  in  the  pressure  traces  taken  by  Bertrand  (13)  or 
Merritt  (45).  It  appears,  therefore,  that  the  ratio  l/l  is  in  correlation 
with  this  pressure  difference;  with  the  increases  in  the  wedge  angle,  the 
bow  wave  becomes  more  and  more  normal  to  the  induced  flow  behind  the 
incident  shock,  so  that  the  pressure  in  region  4  becomes  higher  and  the 
secondary  triple  point  T*  is  "pushed"  upstream.  And  the  ratio  1/ 1,  decreases 
as  a  result. 
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It  seems  possible  that  this  pressure  difference  between  region 
2  and  4  is  a  significant  factor  in  the  genesis  of  the  secondary  Mach  shock 
M' ,  for  when  the  pressure  difference  increases,  some  drastic  mechanism  such 
as  a  shock  is  necessary  to  raise  the  pressure  of  region  2  to  that  of  region 

4.  We  note  that  in  the  DMR  configuration,  the  gas  in  region  5  that  has 
passed  the  straight  reflected  shock  R  is  separated  from  the  gas  in  region  4 
that  has  passed  the  bow  wave  by  the  secondary  slipstream  S'  (see  Fig.  18a). 

The  extent  to  which  the  failure  of  the  assumption  regarding  the 
velocity  of  the  kink  K  (or  the  secondary  triple  point  T')  relative  to  the 
(primary)  triple  point  T  will  have  on  the  transition  boundary  predicted  by 
Ben-Dor's  criterion  requires  further  analysis. 

5.  EXPERIMENTS 

The  experiments  were  performed  in  the  10  cm  x  10  cm  hypervelocity 
shock  tube  and  its  23-cm  dia  plate  Mach-Zehnder  interferometer  at  UTIAS. 

As  the  detailed  accounts  of  the  structures  and  operation  of  the  facilities 
have  been  described  elsewhere  [Boyer  (16), Hall  (28),  and  Ben-Dor  §  Whitten  (10)], 
only  the  special  considerations  relevant  to  the  present  study  will  be  included 
here. 

The  shock  tube  consists  of  a  driver  section  1.8-m  long  and  a 
channel  12-m  long,  with  10  cm  x  18  cm  cross-section,  separated  by  a  diaphragm. 
Mylar  polyester  films  of  various  thickness  were  combined  and  employed  as 
diaphragms.  The  initial  pressure  ratios  across  the  diaphragms  were  computed 
using  the  shock-tube  relations  with  an  adequate  allowance  to  compensate  for 
shock -wave  attenuation  effects  [Boyer  (16)].  Initially,  the  channel  was 
evacuated  by  three  types  of  vacuum  pumps  in  series--mechanical ,  Roots,  and 
diffusion  pumps--down  to  5  to  15  yHg.  The  vacuum  was  measured  by  a  Pirani 
vacuum  gauge  (Type  GP  140).  The  initial  pressures  ranged  from  3.5  torr  for 
high  shock  Mach  numbers  to  80  torr  for  low  shock  Mach  numbers.  Relatively 
high  initial  pressures  were  required  to  meet  the  triggering  threshold  of 
the  pressure  transducers,  which  were  used  for  measuring  the  traversing 
velocity  of  the  shock,  as  well  as  to  maintain  a  reasonable  degree  of  reprod¬ 
ucibility  of  shocks.  Owing  to  air  leaks  into  the  tube  and  the  accuracy  of 
the  vacuum  gauge,  the  purity  of  the  gas  ahead  of  the  shock  was  expected  to 
be  0.04%  to  0.4%,  The  commercially  available  carbon  dioxide  of  purity 
exceeding  99.8%  was  used  as  the  test  gas.  Carbon  dioxide  for  weak  shocks, 
and  helium  or  hydrogen  for  strong  shocks,  were  employed  as  driver  gases. 

The  initial  pressure  was  measured  by  a  U-tube  manometer,  partially  filled 
with  an  oil  of  known  specific  weight,  for  up  to  28  torr  and  by  a  Wallace 
8  Tierman  pressure  gauge  (Model  FA  160-PP10936)  for  over  28  torr.  Possible  errors 
in  the  initial  pressure  do  not  exceed  0.8%.  The  whole  procedure  of  firing 
a  run  after  the  isolation  of  the  channel  took  10  to  25  minutes,  depending  on 
the  pressure  ratio  across  the  diaphragm.  For  lower  Mach  numbers,  the 
reproducibility  of  a  shock  was  quite  good,  while  significant  scatter  proved 
unavoidable  for  higher  Mach  numbers. 

The  initial  temperature  of  the  test  gas  was  taken  to  be  the  same 
as  that  of  the  shock  tube  at  the  midsection  which  was  measured  by  an 
ordinarily  mercury  thermometer  embedded  in  the  tube.  Measurements  of  the 
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temperature  variations  at  the  test  section  by  the  mercury  thermometer  placed 
inside  the  tube  indicated  that  a  temperature  equilibrium  inside  the  tube 
was  achieved  within  5  to  6  minutes  after  the  introduction  of  the  test  gas. 

The  shock  speed  was  calculated  from  the  measured  traverse  time  between  the 
two  consecutive  pressure  transducers  (Atlantic  LD-25) ,  0.305m  apart, 
situated  immediately  upstream  of  the  test  section. 

The  wedges  made  of  steel  of  1.8°,  3.8°,  10.3°,  20.3°,  27.1°, 

30. 1°,  40.2°,  50.4  ,  and  60.1°  semi-vertex  angles  (see  Fig.  37)  were  bolted 
firmly  to  the  bottom  wall  of  the  tube  at  the  test  section.  The  sides  of 
wedges  were  flush  with  the  inside  walls  of  the  tube  and  the  window  glass. 

Although  a  number  of  experiments  were  done  by  using  1.8°  and 
3.8°  wedges,  the  results  could  not  be  used  for  analysis  since  it  proved 
impossible  to  pinpoint  the  shock  confluence  point  exactly  on  either 
shadowgraphs  or  interferograms  due  to  the  weak  reflected  shock  (see  l-'ig.  38a). 

In  the  present  experiments,  infinite-fringe  interferograms  were 
taken  to  investigate  qualitatively  the  density  distribution  of  the  flow 
field,  and  to  furnish  a  mean?  for  a  direct  check  on  the  assumptions  under¬ 
lying  the  simple  theories  of  two-shock  and  three-shock  interactions,  such 
as  uniformity  of  the  flow  properties  in  the  angular  regions  between  the  gas- 
dynamic  discontinuities  at  the  confluence  point.  All  four  mirrors  of  the 
Mach-Zehnder  interferometer  were  adjusted  parallel  to  each  other  so  that  a 
single  fringe  occupied  the  entire  field  of  view  under  a  "no-fiow"  condition. 
With  this  arrangement,  any  fringe  contour  observed  in  the  "flow”  picture 
corresponds  to  a  line  of  equal  density,  or  an  isopycnic  (or  isopycnal). 

As  well,  the  superposition  method  of  flow  and  no-flow  inter¬ 
ferograms  to  produce  isopycnics  was  applied  to  some  experiments  (see  Figs. 

38b  and  11a) .  When  the  two  negatives  of  the  finite-fringe  interferograms 
corresponding  to  "flow"  and  "no-flow"  conditions  are  superimposed,  tiiose 
regions  where  the  fringes  are  just  180°  out  of  phase  become  visible  as  gray 
strips,  or  isopycnics  in  two-dimensional  flow  corresponding  to  those  observed 
in  the  infinite-fringe  interferograms.  However,  as  evident  in  Figs.  38b  and 
11a,  it  is  extremely  difficult  in  practice  to  obtain  high  enough  contrast  to 
make  the  gray  strips  clearly  visible  in  the  superposition  pictures. 

It  is  desirable  to  prepare  tiie  following  four  kinds  of  inter¬ 
ferograms  of  the  identical  experiments  before  drawing  isopycnics  of  a  flow 
field:  (1)  a  white-light  "flow"  inteiterogram;  (2)  a  "no-flow"  finite- 

fringe  monochromatic  interferogram ,  (3)  a  "flow"  finite-fringe  monochromatic 
interferogram;  and  (4)  a  "flow"  infinite-fringe  interferogram.  (1)  determines 
the  fringe  shifts  across  shocks  from  which  density  at  any  point  in  the  flow 
field  may  be  determined  by  (2)  and  (3),  while  (4)  provides  an  aid  in  connect¬ 
ing  the  points  of  equal  density.  Thu?,  the  process  demands  a  high  degree  of 
reproducibility  of  shocks,  since  these  Four  interferograms  should  ideally  be 
taken  under  the  identical  conditions. 

For  small  variations  of  density,  the  infinite-fringe  interferogram 
sometimes  proves  to  be  not  sensitive  enough  to  determine  the  accurate  order 
of  a  fringe  [Whitteman  (66)]. 
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The  adjustments  on  the  interferometer  for  the  infinite-fringe 
interferograms  were  not  without  difficulty:  as  the  four  mirrors  were 
brought  to  a  near  perfect  parallel  to  each  other,  an  "island"  began  to 
appear  in  the  field  owing  to  the  minute  imperfections  in  the  optics  — 
which  are  hardly  noticeable  in  the  finite-fringe  interferograms  —  and  if 
this  "island"  was  successfully  eliminated  from  the  field  of  view,  then  two 
to  three  nearly  parallel  fringes  started  to  appear  in  an  extreme  upper 
corner  instead.  Fortunately,  however,  these  residual  fringes  were  situated 
far  removed  from  the  critical  area  in  the  field,  i.e.  the  vicinity  of  the 
wedge  surface,  their  effects  on  the  quality  of  the  infinite-fringe  inter¬ 
ferograms  were  negligible;  and  all  the  infinite-fringe  interferograms  in 
the  present  study  were  taken  under  that  condition. 

The  shadowgraphs  were  made  using  the  same  equipment  but  obstruct¬ 
ing  the  passage  of  lower  arm  of  the  light  path  by  placing  a  cardboard  in 
front  of  a  mirror. 

The  Control  Data  pulsed  laser  (Model  104A)  was  used  as  the 
light  source.  It  provided  a  monochromatic  light  (69438)  of  approximately 
30  nsec  duration;  and  was  triggered  by  an  i  just  able  delay  circuit  to  make 
it  possible  to  photograph  the  shock  wave  in  the  desired  position.  A  band¬ 
pass  of  the  identical  wavelength  was  placed  in  front  of  the  camera.  If 
desired,  it  was  possible  to  make  interferograms  with  the  second  harmonic 
light  (34728)  simultaneously  with  the  fundamental  mode.  With  this  arrange¬ 
ment,  it  is  possible  to  evaluate  various  effects  that  involve  changes  in 
refractive  indeces  of  two  constituents;  for  example,  dissociation  of  diatomic 
gases,  ionization  of  monoatomic  gases,  etc. 

The  interferograms  and  shadowgraphs  were  photographed  by  using 
4"  x  5"  kodak  Royal  X-Pan  sheet  films  and  enlarged  to  1  to  3  times  the 
actual  size  depending  on  how  detailed  analysis  was  required. 

Two  vertical  threads  1  cm  apart  and  one  horizontal  thread  were 
placed  at  one  side  of  the  window  nearly  2  cm  away  from  the  outer  surface 
of  the  glass  as  the  reference  lines  on  the  pictures  (see  Figs.  38a,  15,  and 
17,  for  example). 

The  Reynolds  number  per  cm  for  the  present  experiments  ranges 
from  5(100  for  lower  Mach  numbers  to  4000  for  the  higher  Mach  numbers. 

6.  CO-i  —  RLLAXA'.  i  .(  I'ITRa/IR  AND  EQUATION'  OF  STATE 

When  a  polyatomic  gas,  initially  in  equilibrium,  experiences  a 
sudden  compression  by  the  passage  of  a  shock  wave,  it  requires  a  length  of 
time  to  re-establish  an  equilibrium  at  a  new  energy  level.  This  phenomenon 
is  clearly  observable  in  the  finite-fringe  interferograms  in  Figs.  4a,  5a, 

6a,  9a,  and  13a  in  the  form  of  gradual  increases  in  the  fringe  shift 
immediately  behind  the  incident  shock.  As  well,  the  superimposed  finite- 
fringe  interferogram  in  Fig.  11a  shows  gray  bands  immediately  behind  the 
incident  shock  and  the  Mach  stem;  and  the  infinite-fringe  interferograms  in 
Figs.  10a,  14a,  and  16a  show  isopycnics  in  the  corresponding  locations. 

The  equilibration  processes  between  various  degrees  of  freedom  of  the 
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molecules  were  treated  by  Bethe  ti  Teller  (12),  who  incidentally  first 
pointed  out  the  possibility  of  the  fully  dispersed  shock  wave.  Kantrowitz 
(38)  investigated  the  relaxation  phenomena  from  the  point  of  heat  capacity 
lag  using  carbon  dioxide  as  a  model.  Griffith  et  al  (23)  found  that  in 
CO,  the  valence  vibrations  adjusted  at  least  100  times  more  slowly  than  do 
the  bending  modes  and  suggested  the  possibility  of  separate  relaxation  times 
for  each  vibrational  mode.  Griffith  and  Kenney  (24)  reported  the  computa¬ 
tional  results  on  partial  equilibrium  state  behind  shock  waves  and  the 
experimental  results  that  excitation  of  each  vibrational  mode  and  dissoci¬ 
ation  took  place  at  rates  which  differed  by  at  least  an  order  of  magnitude 
so  that  the  gas  would  go  through  a  series  of  states  of  partial  equilibrium. 
Schwartz  (52)  developed  the  relaxation  equations  which  applied  to  the 
excitation  of  the  lowest  vibrational  modes  of  CO,  (doubly  degenerate),  and 
predicted  direct  excitation  of  the  bending  mode  and  indirect  excitation  of 
the  valance  modes,  which  was  later  confirmed  by  Witteman  (67)  in  the 
temperature  range  440°  -  816°K.  Johannesen  et  al  (36)  found  that  the 
relaxation  frequency  was  not  a  simple  function  of  translational  temperature 
but  depended  strongly  on  the  departure  from  equilibrium.  Times  of 
vibrational  relaxation  and  dissociation  behind  shock  waves  in  various 
diatomic  and  triatomic  gases  were  measured  by  Gaydon  5  Hurle  (21).  Camac 
(18)  investigat'd  t  nc  structure  ot  shock  waves  in  CO,  in  the  Mach  number 
range  from  5  to  25  and  found  that  all  four  vibrational  modes  relaxed  at 
the  same  rate  and,  that,  at  4  mm  per  microsecond  shock  speed,  the  vibrational 
relaxation  became  comparable  with  the  shock  thickness.  He  found  that  the 
COt  vibrational  relaxation  time  x  was  given  by 

t  -  (1/P)exp(36. 5T  ^//^-3.9)microsecond 

where  P  is  pressure  in  atm  and  T  in  °K.  The  experiments  by  Weihs  fj 
Manheimer-Timnat  (62)  showed  dissociation  relaxation  time  for  CO,  to  change 
from  650  microsecond  at  Mg  =  7.5  to  about  400  microsecond  at  Ms  =  10,  for 
an  initial  pressure  of  3mm  Hg.  (These  happen  to  coincide  with  the  upper 
limit  of  the  present  experiments.)  Significantly,  their  values  were  about 
a  half  of  those  given  by  Gaydon  f(  Hurle  (21).  The  difference  in  the  methods 
of  measurement  was  considered  responsible  for  the  discrepancies.  The  relevant 
parameters  of  Gaydon  Hurle  (21)  and  We*hs  i]  Manheimer-Timnat  (62)  are 
reproduced  in  Table  B1  of  Appendix  B. 

Further  discussion  on  the  thennou\Tiamic  states  and  the  proposed 
models  of  equations  of  state  for  CO,  is  included  in  Appendix  B. 

It  should  be  remembered  that,  by  the  relaxation  time  is  meant 
the  time  required  for  the  ueparture  from  equilibrium  to  fall  to  1/e,  or 
approximately  37%  of  the  equilibrium  value,  so  that  under  the  usual 
assumption  that  the  approach  to  the  equilibrium  is  exponential,  it  takes 
a  considerable  more  time  than  those  quoted  in  Table  B1  to  reach  the  final 
equilibrium  where  no  further  changes  in  the  gas  takes  place. 

The  length  of  the  relaxation  zone  in  the  photograph  can  be 
determined  by  multiplying  the  relaxation  time  in  the  laboratory  frame  of 
reference  by  the  shock  speed.  For  CO,,  in  the  experiments  by  Gaydon  5 
Hurle  (21),  the  vibrational  relaxation  length  ranged  from  0.3  to  3  cm 
and  the  dissociation  relaxation  length  was  approximately  25  cm,  for 
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Ms  -  9;  while  Weihs  §  Manheimer-Timnat  (62)  found  that  the  dissociation 
relaxation  length  to  be  12  cm  at  Mg  =  7.5  and  about  8  cm  at  Ms  =  10. 

The  initial  pressure  in  all  cases  cited  were  approximately  3  torr,  which 
compares  with  the  present  experiments. 

As  well,  it  will  be  recalled  that  the  lower  the  shock  Mach 
number,  the  longer  the  relaxation  times,  as  can  be  observed  in  Fig.  B4 
in  Appendix  B.  This  is  a  direct  consequence  of  lower  temperature  rises 
across  the  shocks  of  lower  Mach  numbers,  and  of  the  fact  that  the  logarithm 
of  the  relaxation  time  is  proportional  to  T“^^,  where  T  is  temperature 
[see  eg.  Kantrowitz  (38),  Camac  (18)].  According  to  Fig.  B4,  the  vibra¬ 
tional  relaxation  zone  behind  the  shock  extends  approximately  10  cm  for 
Mg  =5. 


In  the  light  of  these  data,  and  in  the  light  of  the  fact  that 
the  extent  of  the  zones  of  interest  downstream  of  the  incident  shock,  in 
which  the  shock  angles  are  measured,  is  approximately  1  cm  in  most  cases 
in  the  present  experiments,  it  appears  reasonably  adequate  to  assume  the 
perfect  gas  model  with  y  =  1.290  for  CC>2  for  the  present  study,  the 
prin.  i\  !  e  \  ]-  i  -'i  i  -  delineate  domains  of  various  types  of 
reflection  with  engineering  accuracy. 

7.  PREDICTED  RFGl PA'S  OF  VARIOUS  TYPES  OF  SHOCK  REFLECTION  AND 
COMPARISON  WTTH  EXPERIMENTS 


The  measured  parameters  of  each  e;  ■  sriment  in  the  present  study 
are  listed  in  Tables  1  and  2.  Theoretical  transition  boundaries  of  four 
types  of  reflection  (RR,  SMR,  CMR,  and  DMR)  are  plotted  in  the  shock  Mach 
number  -  effective  wedge  angle  plane  (the  Mg  -  0^  plane)  in  Fig.  34*,  using 

the  five  models  of  equation  of  state  discussed  in  Appendix  B.  The  influence 

of  variation  of  y  op.  the  transition  angles  are  made  manifest  in  Figs.  34a 
through  34e.  Figure  34a  is  based  on  the  perfect-gas  model  with  y  =  1.290 
(Model  P) ,  and  shows  an  excellent  agreement  between  theory  and  experiments. 
This  indicates  that  it  is  necessary  and  sufficient  to  consider  the  gas 
perfect  when  calculating  the  transition  angles. 

Transition  boundary  1  indicates  the  von  Neumann  criterion,  and 
separates  the  Mg  -  '  olane  into  two  domains:  in  the  upper  domain  regular 

reflection  occurs,  while  in  the  lovei  domain  some  form  of  Mach  reflection 

occurs.  The  detachment  criterion  is  indicated  by  boundary  2,  which  also 
divides  the  entire  plane  into  two  domains  in  a  like  manner  as  boundary  1. 

The  overlapped  region  between  boundaries  1  and  2  indicates  the  domain  of 
conflict  of  the  two  criteria,  the  von  Neumann  and  detachment,  and  according 
to  the  former,  RR  is  no  longer  possible  in  this  region.  Boundary  3  is 


*trror  bars  are  not  shown  in  Fig.  34,  since  both  and  Ms  can  be  measured 
to  an  order  of 
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drawn  using  the  criterion  proposed  by  Ben-Dor  5  Glass  (7a, b)  that  the 
transition  from  CMR  to  DMR  occurs  when  the  flow  downstream  of  the  reflected 
shock  is  supersonic  with  respect  to  the  kink  of  the  reflected  shock. 

Boundary  4  corresponds  to  the  condition  NH  =  1,  so  that  the  flow  behind 
the  reflected  shock  is  supersonic  relative  to  t lie  triple  point  T  in  the 
domain  above  boundary  4  and  subsonic  below  it.  Thus,  according  to  the 
definition  of  complex  Mach  reflection  (Sec.  5.2.2),  boundary  4  delineates 
the  domain  of  CMR  from  that  of  SMR .  Boundary  5  indicates  the  limiting 
angles  for  the  given  Ms  at  which  flow  behind  the  incident  steady  shock 
(Fig.  2a)  is  just  sonic  (Mj  =  1),  so  that  no  solution  is  possible  by  means 
of  the  two-shock  or  three-shock  theories  in  the  domain  below  boundary  5 
where  the  flow  behind  the  incident  shock  is  subsonic.  This  boundary  is 
sometimes  called  "the  extreme  sonic  angle"  [see  Smith  (56)]. 

To  the  accuracy  of  the  resolution  (0.2mm)  of  interferograms  in 
the  present  experiments,  the  reflections  in  the  overlapped  regions  between 
boundaries  1  and  2  were  observed  to  be  regular  (see  Figs.  4-6).  An 
excellent  agreement  between  the  predicted  regions  for  CMR  and  DMR  and 
experiments  clearly  confirms  the  validity  of  the  criterion  suggested  by 
Ben-Dor  §  Glass  (7a, b) ,  when  using  the  effective  wedge  angle  0^. 

Figure  35a  shows  the  predicted  domains  of  various  types  of 
reflection  in  the  incident  shock  Mach  number  -  actual  wedge  angle  plane 
(the  Ms-6wplane).  The  actual  wedge  angle  is  obtained  by  subtracting  from 
the  effective  wedge  angle  of  Fig.  35  the  trajectory  angle  x  of  the  triple 
point  which  is  calculated  by  using  equation  5  in  Appendix  D.  Figure  40 
shows  good  agreement  between  the  predicted  domains  and  the  present  experiments. 
Note  that  the  boundary  corresponding  to  the  von  Neumann  criterion  is  not 
shown  in  either  Figs.  35a  or  35b,  since  the  results  of  the  present  study 
support  the  detachment  criterion. 

8.  CONCLUSION 


As  a  continuation  of  the  work  by  Ben-Dor  d  Glass  (7a, b),  who 
investigated  the  transition  boundaries  between  various  types  of  shock 
reflection  in  argon  and  nitrogen,  the  present  study  addressed  itself  to 
the  case  of  a  triatomic  gas  using  OOj  as  the  test  gas.  The  transition 
boundary  map  in  the  Ms-bw  plane  (Fig.  35a),  where  Ms  is  the  shock  Mach 
number  and  0W  the  actual  wedge  angle,  for  CO2,  similar  to  those  reported  in 
Refs.  7a, b  for  the  cases  of  monatomic  and  diatomic  gases,  was  obtained  with 
a  view  to  furnish  a  kind  of  guidepost  for  the  design  of  experiments  in 
future  researches,  excellent  agreement  between  the  boundaries  predicted  by 
Ben-Dor  §  Glass  (7a, b)  and  the  present  experiments  demonstrates  the  validity 
of  the  criteria  unequivocabiy  for  perfect  CO2. 

The  infinite-fringe  interferometry  has  proved  an  excellent 
means  for  investigating  the  density  contours  of  the  flow  fields  in  CO2. 

Some  of  the  infinite-fringe  interferograms  taken  in  the  present  experiments, 
those  of  CMR  and  DMR  in  particular,  may  prove  valuable  when  used  for  compari¬ 
son  with  the  density  contours  predicted  by  means  of  the  numerical  simulation 
of  the  phenomena,  toward  which  a  great  deal  of  effort  is  being  made  at 
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various  research  institutions  at  the  present  moment  [see,  eg.  Book  et  ai  (15J]. 

Of  the  various  models  for  equation  of  state  for  CO2,  the  perfect- 
gas  model  with  y  =  1.290  was  found  to  be  the  most  appropriate  one  to  be  used 
in  conjunction  with  the  simple  two-shock  and  three-shock  theories  over  the 
range  of  shock  Mach  number  from  1.8  to  10.2  of  the  present  experiments. 

Each  of  the  present  experiments  was  analyzed  by  means  of  the 
shock  polar  diagram;  the  results  disagree  with  the  hypothesis  made  by 
Henderson  8  Lozzi  (29)  that  an  apparent  regular  reflection  beyond  the  von 
Neumann  criterion  was  actually  a  nascent  double  Mach  reflection  too  small 
to  be  identified  as  such. 

In  view  of  the  fact  that  there  are  now  available  the  shock 
transition  maps  for  monatomic,  diatomic  (Refs.  7a, b),  and  a  triatomic 
gas,  the  logical  next  step  will  be  to  obtain  a  similar  map  for  air, 
which  is  perhaps  of  the  most  practical  interest.  This  research  is  presently 
underway  at  UTIAS. 
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Si,  s2:  SONIC  POINTS 

(0) ,  (1)  AND  (2)  INDICATE  THE  POINTS  CORRESPONDIN' 
TO  REGIONS  0,  1  AND  2,  RESPECTIVELY. 


FIG.  7  MACH  REFLECTION ,  FLOW  VELOCITIES  IN  THE  FRAME  OF  REFERENCE 
ATTACHED  TO  THE  TRIPLE  POINT  T. 

I:  INCIDENT  SHOCK 

R:  REFLECTED  SHOCK 

M:  MACH  STEM 

S :  SLIPSTREAM 

ew:  WEDGE  ANGLE 

X:  TRAJECTORY  ANGLE  OF  THE  TRIPLE  POINT 

=  0m  +  y:  EFFECTIVE  WEDGE  ANGLE 


FIG.  8  SINGLE  MACH  REFLECTION,  (a)  PHYSICAL  PLANE  (SCHEMATIC), 

(b)  e-P  SHOCK  POLAR  DIAGRAM. 

Is  INCIDENT  SHOCK  POLAR 

R:  REFLECTED-SHOCK  POLAR 

Si,  S2:  SONIC  POINTS  OF  I-  AND  R-POLARS,  RESPECTIVELY. 
NUMBERS  IN  PARENTHESES  INDICATE  THE  STATES  IN  THOSE  REGIONS 


P/PO 


(a) 


THEIR 

(b) 


FIG.  12  COMPLEX  MACH  REFLECTION. 

(a)  PHYSICAL  PLANE  (SCHEMATIC) 

(b)  9-P  SHOCK  POLAR  DIAGRAM 

I:  INCIDENT-SHOCK  POLAR 
R:  REFLECTED-SHOCK  POLAR 
S:  SONIC  POINT  (Mi  =  1) 


(c)  MACH  NU  MBERS  OF  T'."E  FLOW  ALONG  THE  REFLECTING  SURFACE. 


P/PO 


FIG.  18  -  CONTINUED 


18  DOUBLE  MACH  REFLECTION. 


(c)  THE  SYSTEM  OF  REFERENCE  FOR  DMR  ADOPTED  BY 
BEN-DOR  &  GLASS  (7a,  b) . 


(a)  PHYSICAL  PLANE  (SCHEMATIC) 

(b)  0-P  SHOCK  POLAR  DIAGRAM 


I;  INCIDENT -SHOCK  POLAR 

R:  REFLECTED-SHOCK  POLAR 

Si,  Sg:  SONIC  POINTS  OF  THE  I-  AND  R-POLARS, 
RESPECTIVELY 

NUMBERS  IN  PARENTHESES  INDICATE  THE  STATES  IN 
THOSE  REGIONS. 
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FIG.  33  STATIONARY  MACH  REFLECTION 

I:  INCIDENT  SHOCK 
R:  REFLECTED  SHOCK 
M:  MACH  STEM 
S:  SLIPSTREAM 
T:  TRIPLE  POINT 
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34  COMPARISON  OF  PREDICTED  REGIONS  QF  VARIOUS  TYPES  OF  REFLECTION 
AND  EXPERIMENT. 

(a)  PERFECT  COg:  y  =  1.290  (M)DEL  P) 

TRANSITION  BOUNDARIES: 

(1)  THE  VON  NEUMANN  CRITERION 

(2)  THE  DETACHMENT  CRITERION 

(3)  FLOW  BEHIND  THE  REFLECTED  SHOCK  IS  SONIC  RELATIVE 
TO  THE  SECONDARY  TRIPLE  POINT 

(4)  FLOW  BEHIND  THE  REFLECTED  SHOCK  IS  SONIC  RELATIVE 
TO  THE  PRIMARY  TRIPLE  POINT  (M2  =  l) 

(5)  THE  EXTREME  SONIC  ANGLE  (Mi  =  l) 


TYPES  OF  REFLECTION:  0  =  RR;  °  =  SMR;  A  =  CMR;  +  =  DMR 


EFF.  WEDGE  ANGLE  (DEG) 

.0.00  10.00  20.00  30.00  MO. 00  50.00  60.00  70.00  80.00  90.00 


MODEL,  Q1 


+ 


00 


3.00 


5.00 


I - 1 - 1 - 1 - 1 

7.00  9.00  11.00 


INCIDENT  SHOCK  MACH  NO. 


FIG.  31*  "  CONTINUED. 

(b)  QUASI-EQUILIBRIUM  COs  (MDDEL  Ql) . 


MODEL  Q2 


INCIDENT  SHOCK  MACH  NO. 


FIG.  34  -  CONTINUED. 

(c)  QUASI-EQUILIBRIUM  C02  (MODEL  02) 
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FIG.  34  -  CONTINUED. 


(e)  EQUILIBRIUM  CO2  (MODEL  E) 
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FIG.  35(a)  REGIONS  OF  VARIOUS  TYPES  OF  REFLECTION  IN  THE  SHOCK 
MACH  NO.  -  ACTUAL  WEDGE  ANGLE  PLANE  (PERFECT  CO2 
7  =  1.290). 


GURE  3' 


TABLE  1  EXPERIMENTAL  RESULTS:  REGULAR  REFLECTION. 
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SYMBOLS:  EXP.  HO.  =  EXPERIMENT  NUMBER 

Vs  (m/s)  =  INCIDENT  SHOCK  VELOCITY 
Qw  =  WEDGE  ANGLE 

Ms  =  INCIDENT  SHOCK  MACH  NUMBER 

Co  =  ANGLE  OF  REFLECTION 


TABLE  2  EXPERIMENTAL  RESULTS:  MACH  REFLECTION. 
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EXPERIMENT  Nl’MBER 

TYPE  OF  REJECTION:  1  =  SMR,  2  =  CMR,  3  =  DMR 
INCIDENT  SHOCK  VELOCITY  (M./S) 

INCIDENT  SHOCK  MACH  NUMBER:  Ms 
P0  (torr) 

To  f°C) 

WEDGE  ANGLE  (°) 

TRAJECTORY  ANGLE  x  Cf  'THE  TRIPLE  POINT  (°) 
EFFECT IYE  WEDGE  ANGLE  (°),  ^  +  x 

ANGLE  OF  REFLECTION  w  (°) 
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TABLE  B1  RELAXATION  TIMES  AND  LENGTHS  OF  C02- 
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T(l)  MEASURED  MAXIMUM  TEMPERATURE 

T(2)  TEMPERATURE  FOR  NO  DISSOCIATION 

T( 3)  EQUILIBRIUM  TEMFERATURE 

t(L)  LABORATORY  TIME 

T(p)  particle  time 

i  RELAXATION  LENGTH 

*  CALCULATED  BY  THE  PRESENT  AUTHOR  ON  THE  BASIS  OF  THE 
SOUND  SPEED  a  =  270.2  M/SEC  AT  T  =  300°K  FOR  COg 

**  ASYMMETRICAL  STRETCHING  MODE  (CHARACTERISTIC  TEMPERATURE 

03  =  3380. 3°K) 


TABLE  Cl  GLADSTONE-DALE  CONSTANTS  FOR  NITROGEN  [REPRODUCED  FROM 
BEN- DOR  &  WHITTEN  ( 10) ] . 


'-v^^oonat  ant  • 

Ka 

Ka 

Z1 

Z2 

gas 

1 

1 

2 

2 

Nitrogen 

0.2376 

0.328 

0  .246 

0.331 

-0.0904 

-0.0805 

+  0 .03? 

+  1.7* 

+  0.03* 

+  2.5* 

1)  Ai  =  6943^ 

2 )  -  3^71. 
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Kffll,  Kai,  Kng,  Ka2>  Z±  AND  Zg  ARE  ALL  IN  Cm3/G. 


NB: 


SUBSCRIPTS : 
m  =  MOLECULAR 


a  =  ATOMIC 

1  =  WAVELENGTH  *1  =  6y4 

2  =  WAVELENGTH  7g  =  3^7 


■"TYPOGRAPHICAL  ERROR;  THE  CORRECT  VALUE  IS  -0.0850  (AUTHOR) 


APPENDIX  A 


ANALYTICAL  CONSIDERATION 


In  this  appendix,  analytical  aspects  of  the  shock  reflection 
problems  are  reviewed.  Inasmuch  as  only  strong  to  very  strong  shocks  are 
dealt  with  in  the  present  experiments,  the  analyses  of  the  problem  involving 
weak  and  very  weak  shocks  are  precluded  from  consideration  here.  The 
theories  and  the  method  of  solution,  by  which  the  predicted  transition 
boundaries  between  various  types  of  reflection  are  obtained  will  now  be 
explained. 

1.  Pseudo-Stationarv  Flow 


We  first  note  an  important  property  of  the  flow  involving  the 
reflection  phenomenon  of  a  planar  shock  at  a  compressive  corner  in  a  shock 
tube.  From  the  theoretical  point  of  view,  the  walls  of  the  shock  tube  and 
the  reflecting  surface  may  be  regarded  to  extend  semi-infinitely ,  so  that 
no  characteristic  length  scale  is  involved  in  the  flow.  In  addition,  if 
the  flow  is  considered  two-dimensional,  the  entire  configuration  is  observed 
to  grow  linearly  with  time*,  and  every  elementary  fluid  particle  in  the 
flow  moves  away  radially  from  the  corner  with  a  constant  velocity  particular 
to  a  particle  [Bleakney  f)  Taub  (14)]. 

It  is  therefore  possible  to  make  the  observed  non-stationary 
phenomenon  independent  of  time  through  the  transformation 

C  -  x/t  ,  n  =  y/t 

where  x  and  y  are  Cartesian  coordinate:;  with  origin  at  the  corner,  and  t 
is  measured  from  the  instant  the  incident  shock  has  struck  the  corner. 

According  to  the  principle  of  Galilean  invariance,  the  basic 
laws  of  physics  are  identical  in  all  reference  systems  that  move  with 
uniform  velocity  with  respect  to  one  another.  Therefore,  the  two  space 
coordinates  x  and  y  may  be  measured  in  any  frame  of  reference  that  is 
moving  with  a  constant  velocity  relative  to  the  corner. 

The  streamlines  in  the  £  -  n  plane  in  the  cases  of  RR  and  MR 
are  shown  in  Fig.  Al,  which  is  reproduced  from  Jones  et  al  (7).  P  =  (u^.o), 
where  Uj  is  the  induced  flow  velocity  behind  the  incident  shock  relative  to 
the  shock  tube,  indicates  the  point  to  which  streamlines  in  the  region 
behind  the  incident  shock,  undisturbed  by  the  reflected  shock,  are  directed 
in  this  frame  of  reference.  It  should  be  borne  in  mind  that  relative 
orientation  of  each  streamline  to  the  shocks  varies  if  the  frame  of 
reference  is  changed  via  a  Galilean  transformation. 

2.  The  Simple  Two-Shock  and  Three-Shock  Theories 

If  the  frame  of  reference  is  attached  to  the  confluence  point 
of  the  discontinuities,  the  directions  of  the  flow  in  the  angular  regions 

*for  example ,  Gvozdeva  5  Predvodi te leva  (25)  determined  tiie  configuration 
to  be  self-similar  to  an  accuracy  of  up  to  3°o  using  a  high-speed  camera. 
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delineated  by  the  discontinuities  become  uniform  in  the  sense  that  all  the 
streamlines  in  a  region  have  the  same  direction;  and  it  becomes  possible 
to  apply  the  oblique  shock  relations  across  each  shock  subject  to  the 
following  simplifying  assumptions: 

1)  All  the  discontinuities  are  plane  in  the  vicinity  of  the 
confluence  point. 

2)  The  gas  is  ideal  so  that  transport  phenomena  such  as 
viscosity  and  heat  transfer  may  be  neglected. 

3)  The  flow  variables  are  constant  in  each  of  the  angular 
region  around  the  confluence  point. 

The  flow  properties  of  each  region  are  found  algebraically  by 
successive  applications  of  the  oblique  shock  relations  to  each  shock  until 
the  following  boundary  conditions  are  satisfied: 

For  two-shock  interactions,  it  is  required  that  the  flow 
deflections  across  the  incident  and  reflected  shocks  be  equal  and  opposite; 

For  three-shock  interactions,  t  is  lequired 

1)  That  the  total  flow  deflection  across  the  incident  and 
reflected  shocks  be  equal  to  the  flow  deflection  across 
the  Mach  stem; 

2)  That  pressures  across  the  slipstream  be  equal. 

3.  Method  of  Solution 


Normal  shock  relations  [Liepmann  §  Roshko  (44)] 

P  U  =  Pi_U, 

Ha  a  Fb  b 
pa  +  paua  =  pb  +  pbu£ 
ha  +  (l/2)ua  =  hb  +  (1/2) u^ 

where  subscripts  a  and  b  refer  to  the  states  ahead  of  and  behind  the  shock, 
respectively,  u  is  the  component  of  flow  velocity  normal  to  the  shock,  and 
h  enthalpy.  The  equations  of  state  applied  to  each  shock  are  in  the  forms: 

p  =  p(P,T),  h  =  h(P,T) 

where  P  is  the  pressure  and  T  the  temperature. 

In  the  case  of  a  perfect  gas  with  a  constant  specific  heat  ratio,  the  above 
relations  may  be  reduced  to  closed  forms,  the  Rankin-Hugoniot  relations. 

In  the  case  of  real  gases,  however,  recourse  to  the  method  of  successive 
approximation  must  be  taken.  In  the  present  study,  the  iterative  procedure 
in  terms  of  p^  is  continued  until  successive  values  of  p^  are  within  0.05% 
in  the  case  of  the  quasi -equilibrium  models  for  equations  of  state,  and 
0.5%  in  the  case  of  the  complete  equilibrium  model.  Figure  A2  shows  the 
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wave  angle-flow  deflection  angle  diagrams  for  the  various  models  of  equations 
of  state. 


The  real-gas  effect  is  immediately  seen  in  terms  of  the  greater 
flow  deflection  angles  9^  for  real-gas  models  (Models  E,  Ql,  Q2  and  Q3)  in  Fig.  A2. 
For  example,  the  difference  in  between  the  two  extremes  -  the  equilibrium 
model  (Model  E)  for  which  9^  is  greatest  and  the  perfect-gas  model  (Model  P) 
for  which  9^  is  least  for  a  given  incident  flow  Mach  number  Mq  -  amounts 
to  as  much  as  approximately  9°  for  Mo  =  10.  This  is  a  direct  consequence 
of  the  absorption  of  heat  by  the  internal  degrees  of  freedom  of  the  gas  in 
the  real-gas  model.  An  increase  in  the  number  of  degrees  will  result  in 
a  decrease  in  the  temperature  ratio  T^/Tq,  and  hence,  an  increase  in  the 
compression  ratio  P^/Pq  across  the  shock  for  identical  initial  conditions. 

From  the  expression  for  the  flow  deflection  angle  9^  obtainable  from  the 
oblique-shock  relations  after  a  little  manipulation 

0X  =  -  arctan[(pQ/p1)tan0o] 

where  <0j  is  the  flow  angle,  it  is  obvious  that,  for  an  identical  initial 
condition,  the  greater  the  compression  ratio  (i.e.,  the  more  internal  degrees 
of  freedom  are  taken  into  account) ,  the  greater  the  flow  deflection  angle. 

This  effect  may  be  visualized  from  the  standpoint  of  the  flow  velocity,  since 
its  component  parallel  to  the  shock  is  identical  on  either  side  of  the  shock 
irrespective  of  the  model  of  equation  of  state  for  an  identical  initial  con¬ 
dition,  while  its  component  normal  to  the  shock  decreases  according  to  the 
number  of  the  internal  degrees  of  freedom  is  increased. 

These  solutions  are  again  iterated  in  terms  of  0i  in  the  case  of 
two-shock  interactions,  or  in  terms  of  0j  and  03  in  the  case  of  three-shock 
interactions  until 

|l  -  j  <  0.2  x  10'2 

in  the  case  of  two-shock  interactions,  or 

|1  -  (0!  -  02)/e3  I  <  0.2  X  1(T2 

and  -j 

I1  -  P3/P  j  <  0.5  x  10“J 

in  the  case  of  three-shock  interactions,  where 0i,  0^  and  P^  are  the  wave 
angle,  flow  deflection  angle,  and  pressure  of  region  i  (see  Figs.  5  and  10). 
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(a)  EQUILIEKjEIJM  CO2  (MODEL  E) 


(MODEL  Q3  PO  -  1  5 . 0  (TOR)  T0  =  298.0(K) 


I'T  . .  >  -  CONTINUED. 


(U)  QUASI- EQUILIBRIUM  CO2  (MODEL  Q3) 
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APPENDIX  B 


THERMODYNAMIC  STATES  OF  CO., 


With  the  steady  advancement  in  space  exploration  technology  in 
recent  times,  it  is  no  longer  beyond  the  realm  of  possibility  to  send  space¬ 
craft  to  other  planets,  and  CO2  has  been  gaining  a  new  significance;  for 
example,  studies  on  entry  probes  with  pure  CO2  as  a  model  for  Martian  and 
Venusian  atmospheres  have  been  undertaken  in  the  last  decade  [see  Miller 
5  Wilder  (47)]. 

Consisting  of  linear  triatomic  molecules,  CO2  is  one  of  the 
simplest  of  the  polyatomic  gases.  It  is  in  part  owing  to  this  fact  that  a 
considerable  number  of  research  papers  concerning  the  vibrational  relaxation 
phenomena  of  CO^  have  been  published  in  the  past,  some  of  which  were  cited 
in  Section  6. 

CO2  has  four  vibrational  modes:  doubly-degenerate  bending  modes 
(©2  =  960. 2°K),  the  symmetrical  stretching  mode  (£>1  =  1932. 2°K),  and  the 
asymmetrical  stretching  mode  (63  =  3380. 3° K) ,  where  9  is  the  characteristic 
temperature  related  to  the  characteristic  frequency  v  of  oscillation  of  the 
vibrational  mode  by  e  =  hv/kT,  where  h  is  Planck's  constant,  k  is  Boltzmann's 
constant,  and  T  is  the  translational  temperature.  As  can  be  found  from  the 
harmonic  oscillator  approximation  for  the  vibrational  specific  heat, 

Cvib/R  =  (9/T)2exp(0/T)/[exp(0/T)  -  1]2>  the  contribution  to  the  specific 
heat  at  room  temperature  by  the  doubly-degenerate  bending  modes  amounts  to 
0.9R,  where  R  is  the  gas  constant.  (The  valence  modes  contribute  0.07R.) 

Thus,  as  demonstrated  by  Bazhenova  et  al  (3),  it  is  not  realistic  to  set 
the  frozen  specific-heat  ratio  7  of  CO2  to  1.4,  the  value  obtained  when  the 
vibrational  energy  contributions  are  ignored.  Considering  the  fact  that  the 
isentropic  exponent  7^  of  CO2  at  300°K  is  1.286  (Ref.  46),  7  was  set  to  be 
1.290  in  this  study. 

In  conjunction  with  the  simple  two-shock  and  three-shock  theories, 
the  following  models  for  equations  of  state  were  employed  to  compute  the 
thermodynamic  properties  of  the  gas  in  the  regions  around  the  confluence 
point: 

(1)  Equilibrium  Model  (Model  E) 

The  gas  is  considered  to  be  always  in  equilibrium  including  dissociation 
but  not  electronic  excitation  or  ionization  which  are  negligible  over 
the  shock  Mach  number  range  of  the  present  experiments  (1.8  <  Ms  <  10.2). 
Mole  fractions  of  five  species  (CO2,  CO,  O2,  0,  and  C)  were  computed 
by  the  method  of  White  et  al  (37)  on  the  digital  computer.  The  ex¬ 
pression  for  various  partition  functions  and  concomitant  data  of  Allison 
(38),  and  Allison  5  Newman  (39)  were  used.  The  computed  values  were 
checked  against  and  found  to  be  within  a  fraction  of  a  percent  of  the 
corresponding  values  listed  in  Miller  §  Wilder  (40). 

In  this  model,  therefore,  the  gas  is  considered  to  readjust  to  a  new 
equilibrium  state  infinitely  fast  behind  each  shock  around  the  triple 
point. 
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(2)  Quasi-ilqui  librium  Models  (Models  Ql,  Q2,  and  Q3) 

As  seen  in  Table  Bl,  near  the  upper  limit  of  the  shock-Mach-number 
range  of  the  present  experiments,  where  the  effects  of  dissociation 
may  become  significant,  the  dissociation  relaxation  time  is  at  least 
an  order  of  magnitude  larger  than  the  vibrational  relaxation  time. 

It  may  therefore  be  assumed  that  the  vibrational  relaxation  is  complete 
before  dissociation  becomes  significant.  (Dissociation  of  CO2  is 
negligible  below  1900°K.)  Thus,  it  is  meaningful  to  define  quasi¬ 
equilibrium  states  in  which  no  dissociation  occurs  but  one  or  more 
vibrational  modes  are  in  equilibrium  [see,  e.g.,  Griffith  $  Kenney 
(24)].  In  light  of  the  fact  that  the  characteristic  temperature  for 
the  asymmetrical  stretching  mode  is  almost  twice  the  value  for  the 
symmetrical  stretching  mode  which  in  turn  is  double  the  value  for  the 
bending  modes,  three  models  with  various  excited  vibrational  modes  are 
employed,  and  respective  enthalpies  h  for  given  pressure  P  and  tempera¬ 
ture  T  are  obtained  as  follows: 


(a)  Model  Ql  (all  four  vibrational  modes  are  fully  excited). 


h  =  y  RT  +  RT 


202/T  G1/T  03/T 

exp(  e^/T)  -  1  exp(  f^/T)  -  1  exp(03/T)  -  1  _ 


(b)  Model  Q2  (the  asymmetrical  stretching  mode,  which  has  the  highest 
characteristic  temperature,  is  unexcited,  while  the  symmetrical 
stretching  mode  and  the  bending  modes  are  fully  excited). 


h  =  j  RT  +  RT 


2  02/t  e^/T 

exp(02/T)  - 1+  exp(0 j/T)  - 1 


(c)  Model  Q3  (the  valence  modes  are  unexcited,  and  only  the  bending 
modes,  which  have  the  lowest  characteristic  temperature,  are 
fully  excited) . 


h  =  j  RT  +  RT 


292/T 

exp  (  9,/T  -  1 


For  these  models,  the  equation  of  state  for  a  thermally  peifect  gas. 


P  =  pRT 

is  applicable  because  dissociation  is  assumed  absent. 
(3)  Perfect-Gas  Model  (Model  P) 


In  this  model,  the  relaxation  is  considered  to  take  place  infinitely 
slowly,  so  that  the  thermodynamic  states  of  the  gas  in  all  the  regions 
around  the  confluence  point  are  given  by  the  perfect-gas  equations  of 
state  with  a  constant  specific  heat  ratio  7  =  1,290: 
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P  =  pRT 


h  =  tr/Cr  -  i)]RT 

Figures  B-la  through  B-13  show  the  comparison  of  the  predicted 
angle  y  of  the  triple-point  trajectory  by  Eq.  (5)  in  Appendix  D,  using  the 
five  models  of  the  equation  of  state  for  CO2  with  the  measured  angles  %* , 
in  terms  of  the  deviation  Ay  =  %  -  %* .  The  centre  line  indicates  Ax  =  0, 
so  that  the  closer  the  data  points  to  the  centre  line,  the  better  the 
agreement  between  theory  and  experiments.  Our  choice  of  perfect-gas  model 
in  Section  3  seems  to  be  supported  by  Fig.  D-la,  which  shows  the  least 
overall  deviations  of  the  data  points  from  the  centre  line  among  the  five 
models.  Deviations  are  seen  to  be  greater  for  the  models  with  more  internal 
degrees  of  freedom. 

Examination  of  experimental  data  by  Gaydon  5  Hurle  (21)  reproduced 
in  Table  B1  shows  that  the  dissociation  relaxation  time  is  an  order  of  magni¬ 
tude  longer  than  the  vibrational  relaxation  time,  and  that  the  measured  peak 
temperatures  closely  agree  with  the  theoretical  temperatures  calculated  by 
Gaydon  ^  Hurle.  It  seems  therefore  possible  to  preclude  the  equilibrium 
model  from  consideration  for  the  range  of  shock  Mach  number  of  the  present 
studies. 


Figure  B-2a  shows  the  effect  of  increase  of  heat  capacity  due 
to  excitation  of  three  vibrational  modes  and  dissociation.  The  latter 
phenomenon  is  clearly  seen  to  be  pressure-dependent.  Since  even  a  small 
degree  of  dissociation  necessitates  a  large  quantity  of  heat,  temperature  is 
considerably  lower  in  dissociating  gases  than  in  a  perfect  gas.  At  Ms  =  12, 
the  equilibrium  temperature  -  including  dissociation  -  is  less  than  half  of 
the  perfect-gas  temperature.  As  a  consequence  of  the  lowering  of  temperature, 
density  increases  for  real  gases  as  evident  in  Fig.  B-2b.  Density  is  also 
influenced  by  the  initial  pressure  if  dissociation  is  taken  into  account. 

Even  at  Ms  =  12,  the  equilibrium  density  with  Pq  =  1  torr  is  more  than  twice 
as  much  as  the  limiting  value  of  7.9  for  perfect  CO2  with  7  =  1.29;  and  as 
Ms  increases,  density  continues  to  increase  in  real  gases.  However,  as  far 
as  pressure  is  concerned,  behaviour  of  real  gases  does  not  deviate  appreciably 
from  that  of  a  perfect  gas  as  illustrated  in  Fig.  B-2c.  Figure  B-2d  shows 
the  induced  flow  velocity  behind  the  shock,  nondimensionalized  by  the  sound 
speed  in  the  region  in  front  of  the  shock.  The  behaviour  of  the  curves  may 
be  understood  in  relation  to  Fig.  B-2b:  as  density  behind  the  shock  is 
greater  in  real  gases,  the  equation  of  continuity  of  mass  implies  that  the 
flow  velocity  will  accordingly  be  smaller  in  real  gases.  Figure  B-3  shows 
the  theoretical  curves  for  the  ratio  of  Kmole  of  dissociated  CO2  to  Kmole 
of  undissociated  CO2,  or  compressibility  Z.  It  can  be  seen  from  Fig.  B-3 
that  the  effect  of  dissociation  does  not  become  significant  until  Mg  exceeds 
6. 


Figure  B-4  shows  shock  transit  time,  vibration  and  dissociation 
relaxation  times  in  CO2,  which  is  reproduced  from  Fig.  5  of  Camac  (18).  The 
vibration  relaxation  time  of  CO2  at  one  atmosphere  pressure  is  plotted  in 
terms  of  particle  time  against  temperature  in  Fig.  B-5  by  using  the  formula 
of  Ref.  18  (see  Section  6). 
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FIG.  B3  KMDLE  OF  DISSOCIATED  CARBON  DIOXIDE  PER  KMDLE  OF 
UNDISSOCIATED  CARBON  DIOXIDE  (COMPRESSIBILITY,  Z) , 
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FIG.  B4  COg  SHOCK  TRANSIT  TIME,  VIBRATION  AND  DISSOCIATION 
RELAXATION  TIMES  [AFTER  CAMAC  (47)]. 


FIG,  B5  CO2  VIBRATION  RELAXATION  TIME.  IN  ATMOSPHERES 
PRESSURE  -  MECROSECOND  PARTICLE  TIME  IS  PLOTTED 
AS  A  FUNCTION  OF  T-l/3  (°K-l/3)  (AFTER  REF.  47). 


APPENDIX  C 


COMPARISON  OF  DENSITY  CONTOURS 
OBTAINED  BY  DIRECT  AND  INDIRECT  METHODS 


As  mentioned  in  Section  2,  infinite-fringe  interferograms  visualize 
the  lines  of  constant  density,  or  isopycnics,  and  are  perhaps  the  simplest 
direct  means  to  study  the  characteristics  of  density  fields.  When  the  flow 
field  is  not  too  complicated,  it  is  even  possible  to  determine  the  density  of 
each  isopycnic,  for  we  know  the  density  difference  Ap  between  any  two  consecu¬ 
tive  isopycnics  is  related  to  the  wavelength  Xof  the  light  source,  and  the 
So-called  Gladstone-Dale  constant  usually  denoted  by  K,  which  related  the 
density  change  to  the  change  in  refractive  index  of  the  gas,  by 


A p  = 


X 

KL 


where  L  is  the  depth  of  the  test  section  of  the  shock  tube.  Using  the  values 

for  K  listed  in  Table  Cl,  one  can  calculate  Ap  for  the  case  of  N 2 ;  it  is  euual 

to  0.2876  x  10~4  g/cm^  for  X  =  694 08,  and  0.1389  x  10“^  g/cm3  for  X=  3472A. 

In  view  of  these  utilities  of  the  infinite-fringe  interferometry, 
it  was  felt  desirable  to  repeat  the  experiments  reported  in  Ben-Dor  5  Glass 
(5)  taking  the  infinite-fringe  interferograms  to  compare  with  the  isopycnics 
of  the  flow  fields  they  obtained  by  means  of  the  fringe-shift  analysis  of 
finite- fringe  interferograms.  Because  only  two  types  of  reflection  -  RR  and 
SMR  -  are  dealt  with  in  Ref.  5,  further  two  cases  representative  of  CMR  and 

DMR  reported  in  Ref.  7a  were  also  repeated,  and  similar  comparisons  were  made. 

The  results  are  presented  in  this  appendix. 

Figure  C-la  is  the  infinite-fringe  interferogram  corresponding 
to  Case  1  of  Ref.  5.  The  density  contours  obtained  from  fringe-shift  analysis 
by  Ben-Dor  and  Glass  is  reproduced  in  Fig.  C-lc.  "Smearing"  of  isopycnic 
marked  d  in  Fig.  C-la  indicates  the  disturbance  in  the  flow  field  due  to  the 
change  of  flow  direction  behind  the  curved  reflected  shock,  i.e.  the  bow 
wave,  from  that  along  the  shock-tube  bottom  surface  to  the  direction  along 
the  wedge  surface.  The  shapes  of  isopycnics  are  observed  to  be  simple 
immediately  behind  the  bow  shock  where  the  shock  is  almost  normal  to  the 
induced  flow,  and  behind  the  straight  segment  of  the  reflected  shock. 

Figure  C-2a  is  the  infinite-fringe  interferogram  corresponding  to  Case  2  of 
Ref.  5.  Figure  C-2b  is  the  result  obtained  by  Ben-Dor  and  Glass  from  fringe- 
shift  analysis  and  Fig.  C-2c  is  the  numerical  prediction  by  Kutler  8  Shankar 
(41).  The  shape  of  isopycnic  marked  "a"  in  Fig.  C-2a  is  noteworthy  in  that 
neither  Figs.  C-2b  nor  C-2c  managed  to  depict.  Although  the  resemblance  of 
the  shapes  of  the  isopycnics  between  Figs.  C-2a  and  C-2c  is  somewhat  remote, 
the  ranges  of  numerical  values  p/po  assigned  to  isopycnics  in  the  two  figures 
are  in  remarkably  good  agreement,  while  in  Fig.  C-2c  the  highest  density 
region  behind  the  straight  reflected  shock  is  approximately  15%  greater  than 
those  values  in  Figs.  C-2a  and  C-2c.  Figure  C-3a  is  the  infinite-fringe 
interferogram  for  Case  3  in  Ref.  5,  The  density  contours  obtained  from 
fringe-shift  analysis  by  Ben-Dor  and  Glass  and  from  a  numerical  method  by 
Schneyer  (51)  are  reproduced  in  Figs.  C-3b  and  C-3c,  respectively. 
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A  few  words  may  be  said  about  the  number  of  isopycnics  produced 
in  the  density  contours  obtained  from  different  methods.  In  the  infinite- 
fringe  method,  the  difference  in  density  Ap  between  the  two  consecutive 
fringes  is  Ap  =  (KL) ,  so  that  Ap  can  be  made  small  -  and  hence  more  isopycnics 
will  be  produced  -  by  using  the  light  sources  of  shorter  wavelengths.  Com¬ 
pare  Figs.  C-la  and  C-lb,  which  were  taken  simultaneously  in  the  same  experi¬ 
ment.  The  wavelength  of  the  light  source  for  Fig.  C-la  is  one  half  of  the 
wavelength  (the  second  harmonic)  of  the  light  source  for  Fig.  C-lb.  Hence, 
other  quantities  such  as  density  variation,  K,  and  L  being  identical  in  both 
pictures,  there  are  twice  as  many  isopynics  visible  in  Fig.  C-la  as  in  Fig. 

C-lb.  Before  commenting  on  the  number  of  isopycnics  shown  in  the  density 
contours  from  the  finite-fringe  method,  an  outline  of  the  method  used  by 
Ben-Dor  and  Glass  (5,  7a,  b)  is  in  order.  Although  the  principle  of  their 
modus  operandi  is  no  different  from  the  standard  finite-fringe  method  which 
has  long  been  used  in  the  past,  their  innovation  lies  in  Steps  2  and  3 
described  below,  in  which  the  digitizer  was  employed  to  record  the  fringe 
numbers  of  the  points  along  the  fringes  on  the  magnetic  tapes,  and  the  computer 
was  used  to  extrapolate  the  fringe  numbers  at  points  other  than  those  measured 
data  points,  the  processes  customarily  required  a  great  deal  of  tedious  and 
awkward  labour. 

1.  A  point  in  the  "no- flow"  picture  and  the  corresponding  point  in 

the  "flow"  picture  are  chosen  as  origins,  and  the  x-  and  y-axes  are  drawn, 
with  reference  to  which  the  coordinates  of  any  point  in  the  flow  field  may 
be  given. 


2.  Each  fringe  is  "digitized"  by  recording  the  x-  and  y-coordinates 

of  the  points  along  the  fringe  on  the  magnetic  tapes.  Thus,  discrete  points 
in  the  x-y  planes  in  the  "no-flow"  and  "flow"  pictures  are  assigned  particular 
fringe  numbers. 


3.  At  equidistant  points  along  lines  of  constant  x  (or  y)  in  the  x-y 

plane  drawn  in  the  flow  field,  fringe  numbers  Si  in  the  "no-flow"  picture, 
and  S2  in  the  "flow"  picture  are  determined  by  means  of  extrapolation  from 
the  nearest  measured  points.  The  fringe-shift  values  S  at  those  points  are 
then  obtained  by 


S  =  S2  -  si 


from  which  densities  p  are  found  by 


P  = 


P  + 


_A 

KL 


S 


where  q  is  the  initial  density. 


4. 

density. 


Isopycnics  are  drawn  by  connecting  (by  hand)  the  points  of  equal 
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Thus,  it  can  be  seen  that  evaluated  densities  consist  entirely 
of  extrapolated  values.  Since  the  flow  field  is  considered  two-dimensional 
in  tlie  analysis,  extrapolation  is  troublesome  under  the  best  of  circumstances; 
and  the  fact  of  the  matter  is  when  the  flow  field  is  complex,  the  results 
tend  to  be  no  more  reliable  than  sheer  guesses.  The  reliability 
of  the  outcome  from  the  finite-fringe  method,  therefore,  does  not  depend  so 
much  on  how  numerous  the  isopycnics  are  drawn  as  on  how  closely  the  fringes 
are  placed  in  the  original  interferograms .  Moreover,  owing  to  the  expected 
error  in  measuring  fringe  shift  (usually  considered  to  be  0.1  fringes), 
there  is  a  limit  to  how  minutely  the  density  change  measurement  can  be 
achieved  [see,  for  example,  Ben-Dor  §  Whitten  (10)];  therefore,  the  density 
variation  between  two  consecutive  isopycnics  from  the  finite-fringe  method 
should  not  be  less  than  the  limit  thus  imposed,  regardless  of  the  fringe¬ 
spacing  in  the  original  interferograms. 

Infinite-fringe  interferograms  for  CMR  and  DMR,  and  corresponding 
results  from  fringe-shift  analysis  by  Ben-Dor  and  Glass  (7a)  are  shown  in 
Figs.  C-4  and  C-5.  In  either  case,  the  range  of  values  for  p/po  given  in 
Ref.  7a  is  nearly  15%  larger  than  those  determined  from  the  infinite- fringe 
interferograms  of  the  present  experiments.  The  shapes  of  the  isopycnics 
obtainou  using  the  infinite-fringe  method  also  differ  somewhat  from  the 
the  finite- fringe  results. 

The  method  used  by  Ben-Dor  and  Glass  (5,  7a,  b)  enables  the 
densities  to  be  determined  at  discrete  points,  so  that  the  manner  of  connecting 
points  at  which  the  densities  are  identical  to  a  prescribed  tolerance  tends 
inevitably  to  be  speculative,  particularly  when  the  number  of  points  in  the 
flow  field  at  which  density  is  evaluated  is  small.  In  order  to  minimize  the 
inherent  error  in  judgement  and  to  maximize  the  accuracy  of  the  end  result, 
it  seems  essential  to  make  fringe-spacing  of  finite-fringe  interferograms 
as  narrow  as  possible,  and  to  take  as  many  points  on  fringes  as  practicable 
when  digitizing  the  interferograms.  The  reason  is  to  keep  the  interpolation 
errors  small,  since  the  points  at  which  densities  are  to  be  determined  do  not 
necessarily  coincide  with  those  points  at  which  the  fringe  shifts  are  measured, 
i.e.,  the  points  on  fringes,  cumbersome  interpolation  process  becomes  neces¬ 
sary  to  calculate  the  fringe  shifts  at  points  which  do  not  lie  on  any  fringe. 

Despite  its  shortcomings,  the  approach  taken  by  Ben-Dor  and  Glass 
(5,  7a,  b)  must  be  commended  for  being  the  first  to  attempt  to  put  the  digital 
computer  to  good  use  in  order  to  save  a  considerable  amount  of  tedious  manual 
labour  required  in  the  process  of  drawing  isopycnics  by  the  finite- fringe 
method.  It  may  be  regarded  as  a  springboard  for  more  sophisticated  ones 
which  will  undoubtedly  be  developed  in  the  near  future. 
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(c)  ISOPYCNICS  OBTAINED  FROM  FRINGE-SHIFT  ANALYSIS  OF 
FINITE-FRINGE  INTERFEROGRAMS  [AFTER  BEN-DOR  &  GLASS 
(5)].  Ms  =  2.01,  8W  =  63.4%  P0  =  50  torr,  T0  =  291 


FIG.  C2  -  CONTINUED . 

(b)  ISOPYCNICS  OBTAINED  FROM  FRINGE- SHIFT  ANALYSIS  OF  FINITE- 
FRINGE  INTERFEROGRAM3  (A  =  3^72&)  [AFTER  BEN- DOR  &  GLASS 
(6l)].  Ms  =  4.68,  Bw  =  60°,  P0  =  15.31  torr,  T0  =  298. l°K. 

(c)  ISOPYCNICS  FOR  REGULAR  REFLECTION,  NUMERICAL  PREDICTION 
[AFTER  KUTLER  &  SHANKAR  (62)],  02,  7  =  1.4,  Ms  =  4.71, 

Ow  =  60° . 


FIG.  C4  -  CONTINUED . 

(b)  ISOFYCNICS  OBTAINED  FROM  FRINGE-SHIFT  ANALYSIS  OF 
FINITE-FRINGE  INTERFEROGRAMS  [AFTER  BEN-DOR  &  GLASS 
(23a)].  N2,  Ms  =  6. 90,  9„  =  20°,  P0  =  10.12  torr, 

TO  =  295 .8°K. 
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FIG.  C5  -  CONTINUED. 

(b)  ISOPYCMICS  OBTAINED  FROM  FRINGE-SHIFT  ANALYSIS  OF 
FINITE-FRINGE  INTERFEROGRAMS  [AFTER  BEN-DOR  &  GLASS 
(23a)].  N2,  Ms  =  8.06,  0„  =  26.56°,  Po  =  5.11  torr, 

T0  =  298. 2°K. 


APPENDIX  1) 


DERIVATION  OF  THli  EXPRESSION  FOR  Till:  MACH  NUMBER  01-  THE  FLOW 
BEHIND  THE  REFLECTED  SHOCK  RELATIVE  TO  THE  REFLECTING  SURFACE, 
AND  COMPARISON  WITH  EXPERIMENTS 


Referring  to  Pig.  21a,  the  velocity  u2 '  of  the  flow  in  region  2 
relative  to  E,  where  E  is  the  intersection  point  of  the  extension  of  the 
straight  segment  of  the  slipstream  S  with  the  wedge  surface,  can  be  determined 
as  follows:  First,  obtain  the  three-shock  solution  at  the  primary  triple 
point  T  so  that  all  the  flow  properties  of  the  surrounding  regions  are  known. 
Then,  with  the  angled  as  defined  in  Fig.  21a,  the  velocity  w  of  T  relative 
to  E  is  found  by  applying  the  sine  law  to  triangle  OTE  whose  sides  are  propor¬ 
tional  to  the  velocities  of  its  vertices  relative  to  each  other,  and  is  given 
by 


where 


w  =  uf)sinx/sin^  (1) 

f  -  TT  -  (>•  +  05)  (2) 


Now,  U2  and  w  are  in  the  same  direction,  so  that  u2'  is  found  by  a  straight¬ 
forward  application  of  a  Galilean  transformation 


u2  =  u2  -  w  (.3) 


Thus,  the  Mach  number  Mj '  of  the  flow  in  region  2,  relative  to  E  is  given  by 

»  f 

M2  =  u2/a2 

=  u2/a2  -  luy/a0j ia0/a,) (sinx/sin^J  (4) 

=  M2  -  M0(a0/ao.'  (sin%/sinf) 


where  a^  is  the  speed  of  sound  in  region  i,  which  agrees  with  the  expression 
given  by  Gvozdeva  §  Fokeev  (18). 

Equation  (1)  involves  one  parameter  that  is  extraneous  to  the 
three-shock  solution,  i.e.  the  trajectory  angle  x  of  the  triple  point  T. 

As  the  three-shock  solution  pertains  strictly  to  the  local  problem  around  T, 
one  further  assumption  is  required  to  relate  the  shock  angles  at  T  to  the 
"global"  boundary,  i.e.  in  this  case,  the  wedge  surface.  The  most  obvious 
assumption,  and  the  one  we  shall  adopt  here,  is  that  the  extension  of  the 
tangent  to  the  Mach  stem  at  T  meets  the  wedge  surface  perpendicularly.  The 
required  trajectory  angle  x  is  then  simply  given  by 
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X  =  TT/2  -  0  3 


(5) 


where  0,  is  the  flow  angle  incident  on  the  Mach  stem  and  is  given  in  the 
three-snock  solution.  Under  this  assumption  the  wedge  angle  8W  is  related  to 
the  three-shock  resolution  by 


0w  =^3  “*0 


(6) 


where  0o  is  the  angle  between  the  oncoming  flow  uq  and  the  incident  shock  I, 
and  is  a  known  quantity  from  experiment  through  the  relation 


%  =  Vf2  -  (6*  +  jf) 


(7) 


Note  that  the  unprimed  quantities  y  and  9W  in  Eqs.  (5)  and  (6)  indicate  the 
predicted  values,  while  the  primed  ones  y*  and  0W*  in  Eq.  (7)  indicate  the 
measured  values.  Thus,  by  combining  the  above  three  equations  (5),  (6) 
and  (7) ,  it  is  possible  to  see  how  the  two  quantities  Ay  and  A0w  are  related, 
where  Ay  =  y  -  X*  and  A0w  =  "  ew*>  as  follows:  From  Eq.  (5), 

03  =  7T/2  -  y  (F « ) 

Substitution  of  Eqs.  (5')  and  (7)  into  Eq.  (6)  yields 

ew  =  v/2  -  X  -  TT/2  -  0*  -  X*  (8) 

W  W 


This  becomes,  after  rearranging. 


or 


9W  ‘  9W  =  "(X  "  X*)  C9) 

A0W  =  -Ay 


Tlie  transition  boundary  between  CMR  and  DMR  in  Fig.  D1  is  drawn 
according  to  Eq.  (4)  above,  and  all  other  boundaries  are  the  same  as  in  Fig. 

39.  Inadequacy  of  the  transition  criterion  from  CMR  to  DMR  is  apparent, 
confirming  the  remarks  made  by  Bazhenova  et  al  (4)  and  Gvozdeva  5  Fokeev  (27) 
that  the  criterion  of  M2  =  1  relative  to  E  cannot  be  used  to  give  the  transition 
boundary  between  CMR  and  DMR. 
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INCIDENT  SHOCK  MRCH  NO. 

FIG.  D1  COMPARISON  OF  PREDICTED  REGIONS  OF  VARIOUS  TYPES  OF  REFLECTION 
AND  EXPERIMENT  [PERFECT  C02  (7  =  1.29)]. 

TRANSITION  BOUNDARIES: 

(1)  THE  VON  NEUMANN  CRITERION 

(2)  THE  DETACHMENT  CRITERION 

(3)  FLOW  BEHIND  THE  REFLECTED  SHOCK  IS  SONIC 
RELATIVE  TO  THE  WEDGE  SURFACE 

(4)  FLOW  BEHIND  THE  REFLECTED  SHOCK  IS  SONIC 
RELATIVE  TO  THE  PRIMARY  TRIPLE  POINT  (Nfe  =  l) 

(5)  THE  EXTREME  SONIC  ANGLE  (Mj.  =  l) 

TYPES  OF  REFLECTION:  I  =  RR;  o  =  SMR;  A  =  CMR:  +  =  DMR . 


APPENDIX  E 


INFLUENCES  OF  INITIAL  PRESSURES  ON  THE  IMAGES  OF  ISOPYCNICS 
IN  INFINITE-FRINGE  INTERFEROGRAMS 


As  mentioned  in  Appendix  C,  the  density  difference  Ap  between  two 
consecutive  isopycnics  in  the  infinite-fringe  interferogram  is  fixed  for  given 
wavelength  X  of  the  light  source,  Gladstone-Dale  constant  K  of  the  gas,  and 
depth  L  of  the  shock-tube  test  section.  The  density  variation  in  the  flow 
field  is  enhanced  by  the  higher  initial  pressure,  cince  it  will  take  a 
shorter  distance  for  density  to  change  by  Ap  in  a  gas  of  higher  Po  than  of  lower 
Po«  Thus,  by  increasing  Pg,  more  numerous  isopycnics  will  be  produced  in 
the  infinite-fringe  interferograms  for  otherwise  identical  experiments.  It 
can  be  seen  in  Figs.  El  and  E2,  which  show  RR  and  CMR  with  varied  Po,  that 
every  thick,  smeared  isopycnic  in  the  interferograms  with  lower  Po  is  replaced 
by  two  or  more  isopycnics  that  are  narrower  and  distinct  in  the  interferograms 
with  higher  Pq. 

The  vibrational  relaxation  lengths  iv,  the  distance  required  for 
the  departure  from  equilibrium  to  fall  to  1/e,  are  calculated  by  using  Camac's 
empirical  formula  cited  in  Section  6  (p.  16)*  and  shown  in  the  figures. 

The  conclusion  based  on  the  infinite-fringe  interferograms  in  this 
appendix  is  that  it  is  advantageous  to  use  higher  Po  when  investigating  the 
characteristics  of  the  density  field,  for  then  one  can  obtain  a  clearer 
perspective  of  the  field. 

Another  crucial  point  shown  by  these  interferograms  is  that  although 
the  vibrational  relaxation  lengths  vary  by  an  order  of  magnitude,  the  wave 
angles  and  slipstream  angles  do  not  change  for  a  given  Ms  and  Qw.  This  confirms 
that  the  transition  lines  shown  in  Fig.  38a  are  correct  only  for  y  =  1.290. 

The  author  thanks  Dr.  J.  H.  Lee  and  Mr.  R.  Deschambault  for  making 
available  the  infinite-fringe  interferograms  included  in  this  appendix. 


APPENDIX  F 


COMPARISON  OF  THE  PREDICTED  AND  MEASURED  ANGLES  OF  REFLECTION 


The  angle  to  between  the  trajectory  of  the  triple  point  and  the 
reflected  shock,  or  the  angle  of  reflection,  for  each  experiment  in  the 
present  study  was  measured  by  using  a  conventional  protractor.  The  difference 
Aw  between  the  predicted  and  measured  w  for  SMR,  CMR  and  DMR  are  plotted  in 
Figs.  F-la  through  F-le  for  the  five  models  of  equation  of  state  for  CO2  of 
Appendix  B.  Note  in  Fig.  F-la,  Aw  begins  to  increase  beyond  Ms  »  7,  the  range 
of  Ms  over  which  dissociation  becomes  significant  according  to  the  equilibrium 
CO2  model  (cf.  Fig.  B3) ,  while  agreement  between  theory  and  experiments  is 
good  for  quasi-equilibrium  models  Q1  and  Q2  which  neglect  dissociation  but 
take  into  account  the  excitation  of  the  vibrational  modes;  see  Appendix  B 
for  details.  As  seen  in  Fig.  F-le,  the  perfect  gas  model  with  7  =  1.290 
shows  the  least  Aw  among  the  five  models  for  lower  shock  Mach  number  range 
(Ms  ^  3.5),  the  predicted  w  deviates  from  the  measured  angles  significantly 
for  higher  Ms. 

Inspection  of  the  five  diagrams  shows  that  so  far  as  the  parameter 
u  is  concerned.  Models  Q1  or  Q2  appear  to  be  most  consistent  with  experiments. 
It  is,  however,  important  to  note  here  that,  since  it  was  not  possible  to 
accurately  measure  the  desired  angle  at  the  triple  point  owing  to  the  limiting 
resolution  power  of  photographs,  the  measured  angle  is  the  angle  subtended 
by  the  trajectory  of  the  triple  point  and  the  straight  segment  of  the  reflected 
shock,  which  may  or  may  not  be  equal  to  the  angle  in  the  immediate  vicinity  of 
the  triple  point,  say  within  0.2  mm. 
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FIG.  FI  COMPARISON  OF  Ti£E  PREDICTED  AND  MEASURED  ANGLES  OF  REFLECTION. 

(a)  EQUILIBRIUM  COa  (MODEL  E)  **  =  (^predicted  "  ^measured 


HOOEl.  Q| 


FIG.  FI  -  CONTINUED.  QUASI-EQUILIBRIUM  CO2  (MODEL  Ql) 


APPENDIX  G 


BOW  WAVES 


In  the  shock-tube  flow,  a  bow  wave  is  formed  at  a  concave  comer 
even  if  the  flow  behind  the  incident  shock  is  subsonic.  This  occurs  as  a 
result  of  the  reflection  process  of  the  incident  shock  at  the  concave  corner: 
as  the  incident  shock  impinges  on  the  corner,  the  reflected  shock  becomes  a 
bow  wave  and  begins  advancing  into  the  oncoming  induced  flow  behind  the 
incident  shock.  Thus  even  though  the  flow  is  subsonic  relative  to  the  shock 
tube  (this  occurs  when  Ms  =  2.068  for  7  =  1.4  and  Ms  =  1.906  for  j  =  1.29), 
the  bow  wave  advances  with  a  supersonic  speed  relative  to  the  induced  flow 
as  a  result  of  vector  addition  of  the  two  velocities  in  opposite  directions 
and  it  is  only  by  moving  into  the  oncoming  flow  and  so  maintaining  the  super¬ 
sonic  speed  relative  to  it  that  the  bow  wave  can  exist  at  all  in  this  circum¬ 
stance.  Thus,  if  the  induced  flow  is  supersonic  relative  to  the  shock  tube, 
a  bow  wave  may  be  either  attached  or  detached;  and  in  the  latter  case,  it 
will  eventually  assume  an  equilibrium  position  at  a  fixed  distance  away  from 
the  comer.  If,  on  the  other  hand,  the  induced  flow  is  subsonic  relative  to 
the  shock  tube,  a  bow  wave  will  remain  nonstationary  until  it  meets  the  con¬ 
tact  surface  and  refracts.  For  the  supersonic  induced  flow  relative  to  the 
shock  tube,  the  limiting  angles  at  which  a  bow  wave  detaches  from  the  wedge 
can  be  calculated  by  the  methods  available  in  standard  text  books  [see,  for 
example,  Liepmann  §  Roshko  (65)]. 
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APPENDIX  H 


COMPUTER  PROGRAMS 


Included  in  this  appendix  are  the  listing  of  the  computer  programs 
to  calculate  the  transition  angles  of  shock  reflection,  to  draw  the  maps  of 
various  domains  as  shown  in  Figs.  38  and  Dl,  along  with  the  listings  of  the 
programs  to  calculate  the  thermodynamic  properties  of  equilibrium,  quasi¬ 
equilibrium,  and  perfect  CO2,  and  to  draw  P-  shock  polar  diagrams. 

IBM  system  3033  and  the  IBM  system  360-Gould  plotter  combination 
at  the  University  of  Toronto  Computing  Centre  were  used. 

Detailed  analysis  was  done  on  each  experiment  using  the  listed 
computer  programs.  Sample  output  of  two  experiments,  one  for  regular 
reflection  and  the  other  for  complex  Mach  reflection,  are  also  included  in 
this  appendix. 

The  above  may  be  obtained  by  writing  to  Prof.  I.  1.  Glass  at 

UTIAS . 


to  UTIAS,  if  you  require  a  copy.  Available 


Available 


